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Abstract 


This  research  models  and  analyzes  the  inactivation  of  Bacillus  spores  following  a 
radiation  exposure  and  the  process  enacted  by  the  Bacillus  spore  to  repair  the  resulting 
damage.  Irradiation  of  a  spore  and  the  medium  surrounding  the  spore  induces  chemical 
reactions  that  produce  reactive  oxygen  species  (ROS).  This  research  will  consider  the 
reaction-diffusion  of  these  ROS  throughout  the  spore.  These  ROS  can  react  with  the  spore’s 
DNA  and  enzymes  to  degrade  them  to  such  an  extent  that  the  DNA  cannot  be  repaired  or 
replicated,  thus  causing  spore  death.  In  order  to  survive  a  dose  of  radiation,  a  spore  must 
repair  its  damaged  DNA  during  germination.  The  DNA  repair  process  is  dependent  on 
reactions  catalyzed  by  enzymes  that  remain  viable  after  the  radiation  treatment.  Increased 
damage  to  the  enzymes  during  radiation  exposure  effects  the  rate  at  which  the  spore’s  DNA 
is  repaired.  If  the  enzymes  are  damaged  to  such  an  extent  that  they  cannot  complete  the 
DNA  repair  method,  the  spore  will  be  unable  to  reproduce  and  achieve  cellular  outgrowth. 
A  probability  of  survival  model  is  created  based  on  radiation  damage  due  to  the  reaction  of 
ROS  with  the  spore’s  DNA  and  enzymes  and  the  repair  process. 
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MODELING  RADIATION  EFFECTIVENESS  FOR  INACTIVATION  OF  BACILLUS 


SPORES 


I.  Introduction 

The  objective  of  this  research  is  to  model  the  inactivation  of  Bacillus  spores  following 
a  dose  of  radiation.  Irradiation  of  a  spore  and  the  medium  surrounding  the  spore  induces 
chemical  reactions  that  produce  reactive  oxygen  species  (ROS).  These  ROS  can  react  with 
the  spore’s  deoxyribonucleic  acid  (DNA)  and  enzymes  to  degrade  them  to  such  an  extent 
that  the  DNA  cannot  be  repaired  or  replicated,  thus  causing  spore  death.  The  spore’s  DNA 
repair  process  is  dependent  on  reactions  catalyzed  by  enzymes  that  remain  viable  after 
the  radiation  treatment.  If  the  enzymes  are  damaged  to  such  an  extent  that  they  cannot 
complete  the  DNA  repair  method,  the  spore  will  be  unable  to  reproduce  and  achieve  cellular 
outgrowth.  This  research  will  model  radiation  damage  to  the  Bacillus  spore  and  its  resulting 
repair  mechanism  in  order  to  characterize  a  Bacillus  population’s  probability  of  survival. 
This  field  of  study  has  implications  that  range  from  food  and  water  purification  techniques 
to  remediation  of  contamination  from  a  bioweapon  attack. 

Multiple  intentional  and  unintentional  releases  of  Bacillus  anthracis  (B.a.)  into  the 
environment  have  caused  widespread  population  panic.  These  cases  culminated  in  2001, 
following  the  terrorist  attacks  of  11  September,  when  B.a.  spores  were  mailed  through  the 
U.S.  Postal  Service  to  members  of  the  U.S.  Senate  and  to  national  news  agencies.  These 
biological  attacks  resulted  in  22  combined  inhalational  and  cutaneous  cases  of  infection. 
Of  these  22  cases,  five  Americans  died  and  “the  nation  was  terrorized  in  what  became 
the  worst  biological  attacks  in  U.S.  history”  [5].  These  attacks  illustrated  the  prevalent 
need  of  the  U.S.  to  consider  the  detection,  recognition,  environmental  surveillance,  and 
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decontamination  of  B.a.  The  Working  Group  on  Civilian  Biodefense  identified  B.a.  as  a 
biological  agent  that  may  be  used  as  a  weapon  and  “in  worst  case  scenarios,  could  cause 
disease  and  deaths  in  sufficient  numbers  to  gravely  impact  a  city  or  region”  [32].  The  U.S. 
Department  of  Agriculture  considers  B.a.  among  the  fist  of  pathogens  that  could  be  used 
as  a  bioweapon  [29].  In  addition,  the  Center  for  Disease  Control  (CDC)  has  fisted  B.a.  as 
a  Select  Agent  which  has  “the  potential  to  pose  a  severe  threat  to  both  human  and  animal 
health,  to  plant  health,  or  to  animal  and  plant  products”  [87].  And  finally,  the  Department 
of  Defense  considers  B.a.  as  a  threat  to  military  troops  and  describes  it  as  “a  disease  that 
will  kill,  caused  by  a  bacteria  that  already  has  been  used  as  a  weapon  in  America,  and 
that  terrorists  openly  discuss”  [27].  Therefore,  both  political  and  military  factions  within 
the  U.S.  have  stressed  the  need  for  counter  terrorism  weapons  and  biodefensive  systems 
capable  of  destroying  biological  and  chemical  weapons  [107]. 

Weaponized  B.a.  is  not  the  only  cause  for  concern  with  regards  to  the  Bacillus  genus. 
In  August  2003,  five  members  of  the  same  family  were  admitted  to  a  hospital  with  food 
poisoning  symptoms  and  one  child  subsequently  died  due  to  a  food-borne  illness.  Bacillus 
cereus  (B.c.)  was  determined  to  be  the  causative  agent  [25].  The  CDC  estimates  that  one 
in  six  Americans  contract  a  food-borne  illness  every  year  and  3000  of  these  annual  cases 
are  fatal.  Furthermore,  the  most  common  causes  of  these  cases  are  norovirus  and  bacteria 
[32].  Microbial  growth  of  bacteria  is  thought  to  alter  the  character  of  food  during  handling, 
preparation,  transportation,  and  storage.  The  presence  of  bacteria  in  food  products  can  lead 
to  food  spoilage  and  food  poisoning  [2].  B.c.  is  a  food-borne  pathogenic  bacterium  that  is 
finked  to  food  poisoning  and  causes  gastrointestinal  diseases  with  symptoms  ranging  from 
mild  nausea  to  frequent  vomiting.  However,  as  described  above,  contraction  of  the  disease 
produced  by  B.c.  can  be  fatal  and  the  number  of  food  poisoning  cases  caused  by  B.c. 
is  thought  to  be  underreported  [92].  For  both  food  safety  and  counter  terrorism  reasons, 
additional  research  must  be  completed  on  methods  to  inactivate  Bacillus  spores. 
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Further,  recent  events  indicate  that  current  laboratory  processes  to  inactivate  Bacillus 
spores  are  inadequate.  The  Pentagon  claims  that  ’’insufficient  technical  information  in  the 
broader  scientific  community”  was  the  cause  of  live  B.a.  spores  being  disseminated  from 
a  U.S.  Army  laboratory  to  86  labs  around  the  nation  and  seven  foreign  countries  [74].  The 
spores  were  supposed  to  have  been  killed  via  gamma  irradiation  [102]. 

Specifically,  this  research  will  consider  the  inactivation  of  surface-adsorbed  spores  in  a 
humid  environment.  This  scenario  can  be  applied  to  sterilization  efforts  of  equipment  that 
is  contaminated  with  biological  warfare  agents.  For  example,  researchers  are  examining 
military  aircraft  decontamination  methods  in  order  to  prepare  aircraft  for  use  following 
a  biological  attack  [3].  In  addition,  efforts  are  underway  to  investigate  the  microbial 
decontamination  via  UV  radiation  of  spacecraft  surfaces  [56]. 

In  summary,  the  proposed  research  couples  radiation  dose,  damage,  and  repair  in  order 
to  characterize  a  Bacillus  population’s  probability  of  survival.  This  is  done  with  the  intent 
of  modeling  the  inactivation  of  a  population  of  Bacillus  spores  after  it  experiences  radiation 
exposure.  The  specific  proposed  research  objectives  are  listed  below. 

1.1  Research  Objectives 

•  Analytically  model  the  micro-level  damage  within  a  B.a.  spore  due  to  the  reaction- 
diffusion  of  ROS  throughout  the  spore  following  radiation  exposure. 

•  Numerically  model  and  analyze  the  nucleotide  excision  repair  (NER)  mechanism  of 
a  spore’s  damaged  DNA. 

•  Relate  the  damage  and  repair  models  to  determine  a  B.a.  population’s  probability  of 
survival  based  on  the  threshold  killing  mechanisms  of  B.a.  spores. 

1.2  Overview 

The  objectives  previously  defined  will  be  presented  in  the  following  manner.  The 
background  material  necessary  to  model  biological  processes  is  contained  in  Chapter  2. 
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Chapter  3  explains  the  radiation  chemistry  and  production  of  ROS’  that  cause  damage  to 
a  spore’s  DNA  and  enzymes.  It  describes  the  reaction-diffusion  of  these  ROS  throughout 
the  spore  and  presents  a  mathematical  model  of  the  subsequent  damage  to  the  spore’s  DNA 
and  enzymes.  The  development,  numerical  methods,  and  results  of  the  DNA  repair  model 
will  be  presented  in  Chapter  4.  Chapter  5  illustrates  results  of  a  probability  of  survival 
model  for  a  Bacillus  population  based  on  the  relationship  between  radiation  dose,  damage, 
and  repair.  Conclusions,  as  well  as  any  suggested  future  work,  are  contained  in  Chapter  6. 
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II.  Background 


2.1  The  Bacillus  Spore 

This  research  will  discuss  models  that  theorize  the  effects  of  radiation  damage  to 
different  spore-forming  Bacillus  species  such  as  anthracis  and  cereus.  Other  Bacillus 
spore  formers  that  do  not  produce  human  disease  (i.e.  Bacillus  subtilis  (B.s.)  and  Bacillus 
thuringiensis  (B.t.))  are  often  used  as  simulants  for  study  of  spore  inactivation.  Therefore 
it  is  necessary  to  discuss  the  the  validity  of  using  different  species  of  the  Bacillus  genus 
to  characterize  spore  radiation  damage.  All  spores  in  the  Bacillus  genus  share  the  same 
structural  interior  layers.  Bacillus  species  differ  with  regards  to  the  existence  of  the 
exosporium,  which  is  the  outermost  layer  of  some  Bacillus  species.  Specifically,  the 
Bacillus  subtilis  (B.s.)  species  does  not  have  an  exosporium  while  its  closely  related 
species,  Bacillus  anthracis  (B.a.),  Bacillus  cereus  (B.c.),  and  Bacillus  thuringiensis  (B.t.), 
do  [91].  Immediately  inside  the  exosporium  lays  the  spore  coat  and  some  Bacillus  species 
are  believed  to  have  two  spore  coats  instead  of  just  one.  Table  2.1  outlines  the  differences 
between  the  four  Bacillus  species  mentioned  above.  The  absence  or  presence  of  an 
exosporium  is  indicated  by  -  or  -l-  respectively.  Also  included  are  size  parameters  for 
each  species  given  in  units  of  microns  (pm)  [16].  As  discussed  in  the  next  section,  these 


Table  2.1:  The  Bacillus  genus 


Species 

Exosporium 

#  of  Coats 

Diameter  Range  (jum) 

Volume  Range  (pm?) 

B.a. 

-1- 

1 

0.53-1.11 

0.456-0.828 

B.c. 

-1- 

1 

0.76-1.14 

0.543-0.771 

B.s. 

- 

2 

0.41-0.67 

0.127-0.193 

B.t. 

-1- 

1 

0.59-0.96 

0.424-0.574 

5 


structural  differences  do  not  contribute  to  the  spore’s  radiation  resistance.  For  these  reasons 
it  is  assumed  that  radiation  kill  mechanisms  are  similar  across  the  Bacillus  genus. 

2.1.1  The  Structure  of  Bacillus  Spores. 

In  order  to  adequately  model  damage  to  Bacillus  spores  there  must  be  an  understand¬ 
ing  of  the  dormant  spore’s  structure.  From  the  outside  of  the  spore  to  its  center,  the  spore 
layers  are  the  exosporium,  the  spore  coat,  the  outer  membrane,  the  cortex,  the  germ  cell 
wall,  the  inner  membrane,  and  the  core.  These  are  illustrated  in  Figure  2.1  (the  layers  are 
not  drawn  to  scale)  and  each  layer  will  be  discussed  in  detail. 


CELL 
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Figure  2.1:  Structure  of  a  Bacillus  Spore.  Figure  reprinted  with  permission  from  [91]. 
Copyright  2006  John  Wiley  and  Sons. 


The  outermost  layer  of  the  spore,  the  exosporium,  is  made  up  of  proteins  and  the 
function  of  the  exosporium  is  thought  to  be  related  to  the  types  of  cells  that  the  Bacillus 
spore  infects.  The  exosporium  is  readily  permeable  to  small  molecules  and  it  does  not 
contribute  to  the  spore’s  resistance  to  radiation  damage  [2].  The  spore  coat  is  mainly 
composed  of  proteins  but  also  contains  small  amounts  of  carbohydrates  and  lipids.  The 
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coat  acts  as  the  spore’s  first  line  of  defense  against  some  chemical  infiltration  such  as  lytic 
enzymes  that  can  harm  the  cortex  layer.  However,  it  is  not  considered  to  be  a  permeability 
barrier  for  other  small  molecules  such  as  water  molecules,  and  it  does  not  provide  the  spore 
with  any  radiation  resistance  [2,  91]. 

The  next  layer,  the  outer  membrane,  is  an  essential  structure  for  spore  formation 
however  it  does  not  “constitute  a  significant  permeability  barrier  against  small  molecules 
accessing  the  interior  of  the  spore”  [2].  In  addition,  it  does  not  provide  any  spore  resistance 
to  radiation  [91].  Multiple  researchers  have  proved  that  the  exosporium,  spore  coat,  and 
outer  membrane  allow  the  passage  of  water  molecules  [28,  106].  Specifically,  one  of  these 
experiments  by  Westphal  and  others  found  that  an  increase  in  relative  humidity  caused  a 
spore  to  swell  due  to  water  intake. 

Inside  the  outer  membrane  is  the  cortex  which  is  a  peptidoglycan  layer.  Peptidoglycan 
is  a  matrix  of  sugars  and  short  polypeptides  (a  continuous  chain  of  amino  acids  connected 
head-to-tail  through  peptide  bonds)  [99].  The  cortex  functions  mainly  as  a  retaining 
structure  and  reduces  the  water  content  of  the  spore  core  by  allowing  small  molecules 
to  pass  through  it  [48].  The  cortex  is  the  most  hydrated  region  of  the  spore  as  water  makes 
up  48-57  percent  of  the  cortex  [2].  Immediately  interior  to  the  cortex  is  the  germ  cell  wall 
which  is  also  a  peptidoglycan  layer.  The  germ  cell  wall  does  not  contribute  to  the  spore’s 
resistance  to  radiation  damage  [91]. 

The  next  layer  of  the  spore  is  the  inner  membrane  which  is  made  up  of  lipids  and 
proteins.  The  inner  membrane  is  a  strong  permeability  barrier  which  contributes  to  the 
spore’s  resistance  to  deoxyribonucleic  acid  (DNA) -damaging  chemicals.  “Even  water  may 
cross  the  spore’s  inner  membrane  very  slowly”  [91].  Yet  experiments  have  demonstrated 
that  approximately  97  percent  of  core  water  is  transferable  with  its  outside  environment 
and  while  diffusion  of  water  across  the  inner  membrane  is  at  least  two  orders  of  magnitude 
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slower  than  the  membranes  of  vegetative  cells,  the  movement  of  water  through  the  inner 
membrane  does  occur  [96]. 

The  innermost  spore  layer  is  the  core  which  contains  most  of  the  spore’s  enzymes 
and  all  of  its  DNA  as  well  as  other  small  molecules.  Also  found  in  the  core  are  small 
acid-soluble  spore  proteins  (SASPs)  and  dipicolinic  acid  (DPA).  The  SASPs  are  bound  to 
the  spore’s  DNA  and  both  SASPs  and  DPA  play  a  major  role  in  the  spore’s  resistance  to 
ultraviolet  radiation  which  will  be  discussed  in  Section  3.6.2.  DPA  may  have  a  very  minor 
influence  on  spore  resistance  to  ionizing  radiation  by  acting  as  a  scavenger  for  radicals 
produced  by  water  radiolysis  (see  Section  3.3)  [2,  91].  Another  contributing  factor  to 
a  spore’s  resistance  to  radiation  damage  is  its  low  core  water  content.  The  core  water 
constitutes  only  27  to  55  percent  of  the  core  wet  weight  and  the  extremely  low  amount 
of  free  water  greatly  restricts  the  macromolecular  movement  in  the  core.  Specifically, 
the  small  water  content  may  contribute  to  radiation  resistance  by  reducing  the  amount  of 
hydroxyl  radicals  that  are  formed  but  this  fact  has  not  yet  been  proven  [68].  The  effect  of 
the  water  content  in  the  spore  will  be  discussed  further  in  Chapter  3.  The  next  section  will 
concentrate  on  the  structure  of  the  spore’s  DNA. 

2.1.2  The  Structure  of  DNA  [99, 105]. 

The  spore’s  DNA  is  the  genetic  code  of  the  spore  and  transfer  of  its  information  is 
possible  due  to  its  unique  structure.  The  DNA  is  made  up  of  bases  and  two  sugar  phosphate 
backbones  composed  of  a  deoxyribose  sugar  molecules  linked  to  phosphate  groups.  The 
four  bases  are  adenine,  cytosine,  guanine,  and  thymine  and  can  be  visualized  in  Figure  2.2. 
Note  that  thymine  and  cytosine,  the  pyrimidine  bases,  have  single-ring  structures  while 
adenine  and  guanine,  the  purine  bases,  have  double  ring  structures.  In  addition.  Figure  2.2 
shows  the  connection  between  complementary  bases  adenine  and  thymine  via  hydrogen 
bonds  which  are  represented  by  dashed  lines.  The  hydrogen-bonded  adenine  and  thymine 
entity  is  called  a  base  pair.  Similarly,  guanine  and  cytosine  form  a  base  pair.  The  bonding 
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(a)  Adenine  (A)  and  Thymine  (T)  Base  Pairs  (b)  Guanine  (G)  and  Cytosine  (C)  Base  Pairs 

Figure  2.2:  DNA  Base  Pairs.  Figure  reprinted  with  permission  from  [10].  Copyright 
Creative  Commons. 

of  adenine,  a  double-ring  structure,  with  thymine,  a  single  ring  structure,  is  very  important 
to  maintain  the  three  dimensional  structure  of  the  DNA.  The  bonding  of  two  single  ring 
structures  within  a  DNA  strand  would  create  a  divot  in  its  backbone.  Similarly,  if  two 
double  ring  structures  were  paired  together,  there  would  be  a  bulge  in  the  DNA’s  backbone. 
Any  such  conformational  changes  of  the  DNA  structure  must  be  repaired  before  the  DNA 
can  be  accurately  replicated.  Repair  processes  will  be  described  in  Chapter  4.  The  straight 
lines  in  Figure  2.2  represent  covalent  bonds  which  are  stronger  than  hydrogen  bonds.  The 
two  long  backbone  strands  of  the  DNA  wrap  around  each  other  to  form  a  double  helix  (see 
Figure  2.3). 

The  genetic  information  is  encoded  within  the  DNA  by  the  sequence  of  the  bases 
composing  the  DNA  molecule.  The  spore  uses  its  genetic  information  to  synthesis  proteins 
through  a  two  stage  process  of  transcription  and  translation.  Transcription  provides  for 
the  information  from  a  segment  of  DNA  (a  gene)  to  be  communicated  to  ribonucleic  acid 
(RNA).  Translation  is  the  process  by  which  information  from  the  RNA  is  turned  into 
proteins  as  shown  in  Figure  2.4.  During  translation,  three  bases  encode  for  a  letter  in  the 
genetic  language,  called  a  codon.  A  specific  codon  corresponds  to  a  particular  amino  acid 
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Figure  2.3:  The  DNA  Double  Helix.  Figure  reprinted  with  permission  from  [98].  Courtesy: 
National  Human  Genome  Research  Institute 


(labeled  A  site  in  Figure  2.4)  that  is  placed  in  a  unique  location  of  the  protein  polymer 
chain  (labeled  P  site  in  Figure  2.4)  during  its  synthesis. 

Three  bases  of  ordered  information  allow  64  possible  encoding  combinations.  The 
capability  of  the  three-base  letter  is  more  than  sufficient  to  encode  for  20  amino  acids  that 
make  up  the  proteins,  thus  there  exists  redundancy  in  the  code.  As  a  result,  removing  or 
altering  the  third  base  pair  in  the  sequence  of  three  bases  in  a  codon  may  or  may  not  destroy 
the  information  content  identifying  its  corresponding  amino  acid.  However,  it  should  be 
noted  that  removing  the  first  base  more  frequently  causes  problems.  One  incorrect  amino 
acid  can  dramatically  alter  the  function  of  the  protein.  Hence,  damage  to  a  particular  DNA 
base  can  have  a  range  of  effects  depending  on  the  information  content  of  the  codon  or  the 
information  content  of  the  base  within  the  codon  [48].  The  process  by  which  the  dormant 
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Figure  2.4:  The  Translation  Process.  Figure  released  into  the  public  domain  from  [82]. 


spore  is  formed,  how  it  becomes  a  vegetative  cell,  and  the  method  by  which  it  reproduces 
will  be  described  in  the  next  section. 

2.2  Life  Cycle  of  a  Bacterial  Spore 

A  dormant  spore  is  formed  via  a  multi-step  process  called  sporulation  (refer  to  Figure 
2.5).  The  first  step  of  sporulation  is  the  replication  of  DNA  within  the  core.  DNA 
replication  is  an  enzymatically  controlled  process  that  duplicates  the  base  sequence  of  the 
DNA  double  helix.  Replication  begins  with  the  separation  of  the  two  sugar  phosphate 
backbones.  In  order  to  accomplish  this  separation,  the  hydrogen  bonds  between  the 
complementary  base  pairs  are  severed  while  the  backbones  remain  intact.  A  new  DNA 
strand  is  synthesized  from  the  base  pairs  with  the  cleaved  hydrogen  bonds.  DNA  replication 
is  semiconservative  as  the  two  resulting  double  helices  consist  of  one  old  sugar  phosphate 
backbone  and  one  new  complementary  strand  [54]. 
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After  DNA  replication,  sporulation  continues  with  the  formation  of  the  spore  septum 
which  is  an  ingrowth  of  the  inner  membrane.  The  spore  septum  separates  the  replicated 
DNA  and  a  small  portion  of  the  core’s  cytoplasm  to  include  enzymes.  The  structure 
enclosed  by  the  spore  septum  is  called  the  forespore.  The  next  step  of  sporulation  is  the 
creation  of  a  thick  spore  coat  of  proteins  around  the  forespore.  This  new  structure  is  called 
the  endospore  and  once  it  matures,  the  endospore  is  freed  from  the  vegetative  cell  by  the 
rupturing  of  the  spore  cell  wall  [99] . 


Figure  2.5:  Sporulation.  Figure  reprinted  by  permission  from  Macmillan  Publishers  Ltd: 
Nature  Reviews  Microbiology  [57],  copyright  2012. 


A  dormant  spore  can  survive  for  extremely  long  periods  of  time  and  in  harsh 
environmental  conditions.  In  a  nutrient  rich  environment,  the  endospore  will  undergo 
germination  which  is  the  first  step  in  the  process  by  which  bacteria  transforms  from  a 
dormant  spore  into  a  vegetative  cell  [28].  The  nutrient  induced  germination  process  allows 
changes  to  the  permeability  of  the  inner  membrane.  During  the  first  few  minutes  of 
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germination,  the  spore  core  takes  in  water  and  releases  DPA  [20].  The  water  content  of 
growing  cells  constitutes  approximately  80  percent  of  its  wet  weight.  In  order  to  achieve 
this  exchange  of  molecules  during  germination,  the  inner  membrane  expands  (see  Figure 
2.5).  This  contributes  to  an  increase  in  the  spore’s  susceptibility  to  radiation  damage  during 
germination  because  of  the  increased  molecular  movement  and  amount  of  water  molecules 
in  the  core.  In  addition,  this  process  further  demonstrates  that  the  inner  membrane  is 
permeable.  During  germination  the  cortex  and  spore  coat  are  degraded  by  enzymes  which 
allows  hydration  of  the  spore  core.  Mutants  lacking  a  single  enzyme  are  still  able  to 
germinate,  however  dormant  spores  lacking  both  of  the  two  cortex  degrading  enzymes 
are  not  able  to  initiate  outgrowth  [2].  Rehydration  of  the  spore  core  allows  protein  mobility 
which,  in  turn,  permits  enzyme  action.  This  triggers  the  transition  from  germination  to 
vegetative  metabolism  and  thus  cellular  outgrowth. 

It  is  only  during  germination  that  the  spore’s  enzymatic  activity  allows  the  spore  to 
repair  any  damage  incurred  during  dormancy  [88].  DNA  repair  is  conducted  via  spore 
enzymes  and  will  be  discussed  further  in  Chapter  4.  However,  since  dormant  spores  do  not 
synthesis  enzymes,  the  spore  must  contain  enough  viable  enzymes  to  complete  the  repair 
process.  If  too  much  damage  is  accumulated,  “this  damage  can  overwhelm  the  capacity  of 
repair  systems  and  lead  to  the  death  of  the  germinated  spore”  [90]. 

Note  that  the  sporulation  process  involves  one  vegetative  cell  creating  one  dormant 
spore  thus  sporulation  is  not  a  means  of  reproduction.  Reproduction  is  the  development  of 
a  vegetative  cell  instead  of  an  endospore.  Bacteria  reproduce  through  binary  fission  which 
is  the  formation  of  two  daughter  cells  from  a  single  parent  cell.  Specifically,  the  DNA 
replicates  itself  and  the  germ  cell  wall  and  inner  membrane  elongate  and  begin  to  form 
around  each  of  the  two  chromosomes.  Eventually  the  cell  wall  and  inner  membrane  fully 
encapsulate  each  DNA  strand  and  the  cells  separate  as  shown  by  Figure  2.5  [99]. 
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Outgrowth  of  a  colony  of  bacteria  occurs  in  four  phases;  the  lag  phase,  the  exponential 
phase,  the  stationary  phase,  and  the  decline  or  death  phase.  These  phases  are  depicted  in 
Figure  2.6.  During  the  lag  phase,  the  bacteria  population  experiences  a  large  amount  of 
metabolic  activity  in  which  it  is  repairing  DNA  and  replicating  both  DNA  and  enzymes. 
There  exists  little  to  no  binary  fission  during  this  phase.  It  is  during  the  lag  phase  that 
the  population  is  the  most  susceptible  to  radiation  and  chemical  damage  because  adverse 
conditions  will  interfere  with  the  growth  process.  The  exponential  phase  is  the  phase  in 
which  the  cells  begin  to  divide  and  the  population  undergoes  exponential  growth.  When  the 
production  of  new  cells  is  balanced  with  the  number  of  cell  deaths,  the  colony  of  bacteria 
has  entered  the  stationary  phase.  During  this  period  of  equilibrium,  the  metabolic  activity 
of  individual  cells  is  slowed.  The  death  phase  indicates  that  the  number  of  cell  deaths 
exceeds  the  number  of  vegetative  cells  being  formed.  The  population  size  will  decrease  at 
a  logarithmic  rate  until  only  a  fraction  of  its  cells  remain  or  the  population  entirely  dies  out 
[99]. 


Figure  2.6:  Phases  of  Outgrowth.  Figure  reprinted  with  permission  from  [84].  Copyright 
1987,  American  Society  for  Microbiology. 
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2.3  Summary 

This  chapter  has  described  the  structure  and  life  cycle  of  a  Bacillus  spore  and  the 
structure  of  its  DNA.  These  characteristics  will  play  an  important  role  in  the  damage 
mechanisms  and  ultimately  inactivation  of  Bacillus  spores.  This  research  will  model 
damage  to  the  Bacillus  species  by  irradiation  which  is  discussed  in  the  next  chapter. 
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III.  Radiation  Damage 


Radiation  damage  to  biological  matter  can  be  categorized  into  two  broad  groups; 
direct  and  indirect  action.  Direct  action  is  the  direct  deposition  of  energy  into  a  biological 
molecule  of  interest  such  as  DNA.  A  single  event  of  direct  radiation  action  upon  an 
important  biological  molecule  will  most  likely  result  in  a  loss  of  biological  activity  of 
the  molecule  [4].  This  research  will  focus  on  the  inactivation  of  biological  spores  due  to 
indirect  action. 

Indirect  action  is  “the  interaction  of  solute  molecules  and  the  reactive  species  of 
solvent  molecules  formed  by  the  direct  action  of  radiation  on  the  solvent”  [4] .  Most  indirect 
action  in  aqueous  media  will  result  from  interactions  of  important  biomolecules  in  the  spore 
with  reactive  oxygen  species  (ROS)  formed  by  the  irradiation  of  adsorbed  water  molecules. 
Specific  products  formed  by  the  irradiation  of  water  molecules  are  discussed  in  Section 
3.3.  Prior  to  that,  the  two  types  of  radiation  (ultraviolet  (UV)  and  ionizing  radiation)  will 
be  examined  along  with  the  differences  in  their  chemical  reactions  due  to  various  levels  of 
energy  deposition. 

3.1  Ionizing  Radiation 

Ionizing  radiation  produces  highly  energetic,  ionizing  rays  and  radiolysis  is  the 
absorption  of  a  particle  of  ionizing  radiation.  When  a  quanta  of  ionizing  radiation  interacts 
with  matter,  a  fast  electron  is  produced  which  then  proceeds  to  produce  ionization  along  its 
track.  When  a  material  experiences  radiolysis,  one  quantum  or  particle  of  ionizing  radiation 
excites  many  molecules  within  the  medium  and  the  reactive  species  produced  by  ionizing 
radiation  are  not  homogeneously  distributed  in  the  solution. 

There  are  many  forms  of  ionizing  radiation  such  as  /?-,  y-,  and  X-rays  and  electrons, 
protons,  and  or-particles.  Ionizing  radiations  are  often  categorized  according  to  their  linear 
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energy  transfer  (LET)  [66].  LET  for  charged  particles  of  ionizing  radiation  is  the  rate  at 
which  energy,  E,  is  transferred  to  an  absorbing  material  per  unit  length  of  track,  x,  which 
can  be  represented  mathematically  by 


LET  =  - 


d.r  ’ 


At  lower  energies,  the  Bethe  formula  for  stopping  power,  or  energy  loss  suffered  by  a 
charged  particle,  is  given  by  [64] 


dv: 


A-nz^e"^ 

mv^ 


NZ  In 


2mv^ 

I 


where 


ze  =  charge  of  the  incident  particle 
e  =  electron 

V  =  velocity  of  the  incident  particle 
m  =  mass  of  the  electron 
N  =  molecules  per  volume 
Z  =  number  of  electrons  per  molecule 
I  =  average  ionization  potential. 

As  shown  by  this  equation,  LET  of  lower  energy  particles  increases  as  kinetic  energy 
decreases.  In  other  words,  as  the  ionizing  particles  are  slowed  down  in  the  medium, 
their  LET  increases.  Therefore,  towards  the  end  of  a  track,  the  LET  will  be  much  higher 
because  the  energy  of  the  particle  decreases  significantly.  This  means  that  the  region  of 
ionization  at  the  end  of  the  track  will  be  both  bigger  and  denser.  This  region  is  referred  to 
as  a  blob.  In  addition,  spurs  caused  by  ionization  of  secondary  electrons  may  extend  from 
the  primary  ionization  track.  Eigure  3.1  depicts  the  extension  of  spurs  from  the  main  linear 
path  of  an  ionizing  particle.  The  specific  nature  of  each  track  of  ionization  is  dependent 
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Figure  3.1:  Track  and  Spur  Depiction.  Figure  reprinted  from  [40]  Copyright  (1973),  with 
permission  from  Elsevier. 


upon  the  type  of  ionizing  particle.  For  example,  a  1  megaelectronvolt  (1  MeV)  electron, 

a  subatomic  particle  with  a  negative  charge  which  belongs  to  the  low-LET  class,  ionizes 
approximately  one  molecule  out  of  every  5000  along  its  track.  In  comparison,  a  1  MeV 
cr-particle,  a  positively  charged,  high  mass  particle  with  high-EET,  ionizes  every  molecule 
along  its  track.  This  makes  a-radiation  the  most  destructive  form  of  ionizing  radiation  [40] . 

Eet  A  represent  one  molecule  of  matter.  A  typical  chemical  reaction  resulting  from 
radiolysis  is  represented  by 

A  +  hv  ^  A'^  +  e~  +  hV 

where  hv  is  the  energy  of  the  incoming  photon  and  e  is  one  electron.  This  reaction  describes 
the  Compton  effect  in  which  a  photon  collides  with  an  electron  and  drives  the  electron  out 
of  its  molecule.  The  photoelectric  effect  occurs  when  the  photon  collides  with  an  electron 
and  is  absorbed  by  that  electron.  In  this  case,  the  electron  is  known  as  the  photoelectron 
and  is  denoted  by  e~ .  The  reaction  given  by  the  photoelectric  effect  is  characterized  by  [64] 

A  +  hv  ^  A'^  e~ . 


3.2  Ultraviolet  Radiation 

UV  radiation  emits  quanta  of  energies  (or  photons)  anywhere  from  one  to  ten 
eV  although  it  is  more  commonly  categorized  by  its  wavelength  which  ranges  from 
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approximately  100  to  400  nanometers  (nm).  Table  3.1  distinguishes  between  the  UV-A, 
UV-B,  and  UV-C  ranges  of  the  UV  spectrum. 


Table  3.1:  UV  Categories  [2] 


UV  Type 

Wavelength 

UV-A 

315-400  nm 

UV-B 

280-315  nm 

UV-C 

100-280  nm 

UV  radiation  produces  photolysis,  otherwise  known  as  photochemistry,  in  which 
chemical  reactions  are  caused  by  absorption  of  light.  In  photochemistry,  each  quantum 
of  radiation  is  absorbed  by  one  molecule  which  produces  one  excited  state  therefore  the 
excited  states  are  homogeneously  distributed  in  a  homogeneous  solution.  The  typical 
chemical  reaction  resulting  from  photolysis  is 

A  +  hv  ^  A* 

where  A*  represents  the  excited  state  of  the  molecule  of  matter,  A  [40].  The  specific 
interactive  species  produced  by  the  radiolysis  and  photolysis  of  water  are  discussed  in  the 
next  section. 

3.3  Radiation  of  Water 

The  reactive  species  produced  by  the  ionizing  radiation  of  water  are  OH,  H,  H2O2, 
H2,  H30'^,  and  OH”.  The  chemical  reactions  that  produce  these  interactive  species  will 

now  be  explored.  According  to  Moore  [64],  regardless  of  the  source  of  radiation,  the  first 
two  chemical  reactions  due  to  radiation  are 

H20^H20^+e-^  (3.1) 
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H20^  +  H2O  +  OH  (3.2) 

which  generates  the  hydrated  electron,  and  the  hydroxyl  radical,  OH.  Also  according 
to  Moore  [64],  radiolysis  will  produce  many  OH  radicals  close  together  and  therefore  they 
may  combine  to  produce  hydrogen  peroxide: 

lOH  H2O2. 


A  subsequent  chemical  reaction  is 

OH  +  H2O2  HO2  +  H2O 

which  produces  the  peroxy  radical  HO2.  The  hydrated  electron  from  Equation  3.1  reacts 
with  water  to  create  the  hydrogen  radical: 

e-^  +  H2O  H  +  OH-. 

The  hydrogen  radicals  can  combine  with  each  other  in  the  following  reaction: 

2H  H2. 

The  previous  reactions  are  just  a  few  of  the  possible  reactions  that  occur  in  pure  water. 
If  there  exists  a  thin  layer  of  adsorbed  water  on  a  surface,  then  the  adsorbed  water  would 
be  saturated  with  dissolved  02-  Thus  the  chemical  reaction 

H  +  02^  HO2  (3.3) 

could  occur.  Similarly,  the  presence  of  oxygen  can  produce  the  superoxide  radical  as  shown 
by  the  following  reaction: 

^aq  +  02^  O2. 

All  of  these  reactions  are  governed  by  rates  of  reactions  which  will  be  discussed  in  the  next 
section. 
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3.4  Reaction  Kinetics 


The  study  of  rates  of  chemical  processes  is  known  as  reaction  kinetics.  The  main 
factors  influencing  the  reaction  rates  are;  the  concentration  and  physical  state  of  the 
reactants,  the  temporal  behavior  of  the  reactants  and  the  presence  of  any  catalyst  in  the 
reaction  [67].  The  rate  at  which  HO2  is  produced  from  Equation  3.3  is  given  by 


d[HQ2] 

dt 


=  fc[H]  [O2] 


where 


[HO2]  =  concentration  of  HO2  (M) 
t  =  time  (seconds  (s)) 
k  =  reaction  rate  coefficient  (M“^s“^) 
[H]  =  concentration  of  H  (M) 

[O2]  =  concentration  of  O2  (M). 


The  units  of  [HO2],  [H],  and  [O2]  are  given  by  molarity,  or  molar  concentration.  Molarity 
is  denoted  by  M  and  is  defined  as  the  number  of  moles  of  solute  per  litre  of  solvent  (i.e. 
M  =  mol/L).  The  concept  of  reaction  kinetics  will  be  used  throughout  the  development  of 
mathematical  models  for  damage,  repair,  and  kill  probability. 

Table  3.2  contains  the  rate  coefficients  given  in  units  of  M“^s“^  for  the  chemical 
reactions  described  above  and  many  more.  When  considering  conflicting  reaction  rate 
coefficients  found  in  the  literature,  the  rate  corresponding  to  a  reaction  in  the  most  neutral 
pH  environment  was  recorded.  The  reason  for  this  is  because  most  bacteria  survive  an 
environment  with  pH  measurements  between  five  and  eight  and  “microbial  cells  require 
their  internal  pH  to  remain  constant  in  order  to  maintain  essential  cell  functions”  [93].  It 
should  also  be  noted  that  the  chemical  reactions  that  occur  in  highly  alkaline  solutions  were 
not  considered.  And  finally,  the  reaction  H  +  H2O  OH  +  H2  was  not  included  because 
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Table  3.2:  Reactants  and  Rate  Coefncients  for  the  Radiolysis  of  Aqueous  Solutions  [30, 40] 


Reactant 

Rate  Coefficient  (M  ^5 

H20^  +  H2O  +  OH 

1x10^° 

20H  H2O2 

5.3  X  10‘’ 

OH  +  H2O2  HO2  +  H2O 

3.8  X  10^ 

e-^  +  H2O2  ^OH  +  OH- 

1.3  X  10^° 

e-^  +  H^O^  ^H  +  H20 

2.06  X  10^° 

e-g  +  H2O  H  +  OH- 

1.6  X  10^° 

e-^  +  H^  H2  +  OH- 

3  X  lOi** 

^aq  +  ^aq  ^  ^2  +  26>H“ 

5x  10‘^ 

^aq  +  02^  0- 

1.9  X  10^° 

e-^  +  OH  ^  OH- 

3  X  10^° 

(^~aq  +  ^02  HO- 

2x  10^° 

OH  +  H  ^  H2O 

3.2  X  10^° 

OH  +  H2  ^  H  +  H2O 

6x  10^ 

OH  +  HO2  -^02  +  H2O 

7.1  X  10‘’ 

H  +  H  ^  H2 

7.9  X  10‘’ 

H  +  HO2  H2O2 

2x  lOi** 

H  +  02^  HO2 

2.6  X  10^0 

H  +  H2O2  OH  +  H2O 

3.44  X  10’ 

OH +  0-^02  +  OH- 

9.96  X  10‘' 

OH  +  HO2  HO2  +  OH- 

5x  10‘^ 

HO2  +  O2  — ^  O2  +  HO 2 

9.5  X  10’ 

2HO2  — ^  H2O2  "1"  O2 

8.1  X  10^ 

this  reaction  only  occurs  in  the  absence  of  oxygen  [4].  This  table  is  not  all  inclusive  as 
many  other  chemical  reactions  are  possible;  however,  these  are  the  predominant  reactions. 

The  rate  of  reaction  is  limited  by  the  encounter  rate  of  reacting  species  and  the  rate  of 
mutual  diffusion.  The  diffusion  controlled  rate  limits  the  frequency  at  which  molecules  can 
collide  with  each  other  in  a  solution.  Diffusion  controlled  reactions  occur  so  quickly  that 
the  rate  of  reaction  becomes  the  rate  of  transport  of  the  reactants  through  the  solution.  The 
diffusion  controlled  limit  in  water  is  on  the  order  of  10*°  M“*s“*  which  means  that  in  water 
the  rate  of  reactions  from  Table  3.2  cannot  be  larger  than  10*°  M“*s“*  [105].  Molecular 
diffusion  will  be  discussed  in  the  next  section. 

3.5  Molecular  Diffusion 

Molecular  diffusion  is  the  movement  of  molecules  from  a  region  of  higher  concentra¬ 
tion  to  that  of  a  region  with  lower  concentration.  This  transport  of  molecules  is  founded  on 
random  molecular  movement.  Diffusion  is  mathematically  modeled  via  Pick’s  first  law  of 
diffusion  which  states  that  the  flux  of  a  species  is  proportional  to  its  gradient  of  concentra¬ 
tion  [67]: 

J  =  -DVc  (3.4) 

where 

J  = 

D  = 
c(x,t)  = 

X  = 

t  = 

The  diffusion  coefficient  measures  “how  efficiently  the  particles  disperse  from  a  high  to 
a  low  density”  [67].  The  negative  sign  reflects  the  fact  that  molecular  flow  is  from  a 


flux 


M 

m^  s 


diffusion  coefficient 


m" 

s 


concentration  of  a  species  (M) 
position  (meters  (m)) 
time  (s). 
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high  concentration  to  a  low  concentration.  Table  3.3  contains  some  of  the  molecular 
compounds  produced  by  ionizing  radiation  and  their  diffusion  coefficients  in  water.  All 
diffusion  constants  are  recorded  in  units  of  m^s“^  and  with  water  temperatures  between  20 
and  25  degrees  Celsius  (room  temperature). 

Table  3.3:  Diffusion  Coefficients  [31,  40,  95,  105] 


Reactant 

Diffusion  Coefficient  (m^s 

H2O2 

1.4  X  10-'^ 

HO2 

2  X  10-*^ 

3.62  X  10-^ 

OH- 

1.97  X  10-^ 

OH 

2  X  10-*^ 

H 

4.5  X  10-'^ 

O2 

2  X  10-‘^ 

p~ 

^  aq 

4.25  X  10-^ 

3.6  Radiation  Damage  to  a  Spore 

Another  difference  between  ionizing  and  UV  radiation  is  the  type  of  damage  they 
inflict  upon  a  dormant  spore.  The  next  two  sections  will  describe  the  broad  range  of  damage 
to  the  molecules  within  a  spore  caused  by  ionizing  and  UV  radiation. 

3.6.1  Ionizing  Radiation  Damage. 

Some  of  the  ROS’  discussed  in  Section  3.3  cause  indirect  damage  to  the  spore’s 
DNA.  They  can  produce  radiation-induced  base  modifications  within  the  DNA  such  as 
the  severing  of  an  amino  group  from  a  DNA’s  base  pair.  This  type  of  base  modification  is 
called  deamination.  Another  base  modification  example  is  depurination  which  is  the  loss  of 
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a  purine  base.  Base  modifications  will  most  likely  induce  a  base  pair  mismatch  which  will 
create  a  bulge  or  divot  in  the  DNA  backbone.  Without  repair,  this  structural  conformation 
of  the  backbone  will  prohibit  the  DNA  from  being  read  during  the  replication  process  which 
can  be  lethal  to  the  spore.  Table  3.4  gives  the  rate  of  reaction  of  the  hydrated  electron, 
with  the  DNA  bases.  It  has  also  been  shown  that  base  modifications  can  result  from  OH 
reacting  with  DNA  bases  and  H  reacting  with  the  pyrimidine  bases  [14].  For  example,  the 
reaction  of  OH  with  thymine  will  result  in  an  abstraction  of  an  H-atom  [26]. 

Table  3.4:  Rates  of  Reaction  of  e~  with  DNA  Bases  [66] 


DNA  Base 

Reaction  Rate  (M  ^s 

Adenine 

3x  10^° 

Cytosine 

7-10  X  lO*^ 

Guanine 

2.5  X  10* 

Thymine 

1.7x  lO^o 

In  addition  to  base  pair  reactions,  ROS  can  target  the  DNA  backbone  resulting  in 
a  broad  range  of  radical-induced  products.  For  example,  the  OH  radical  can  react  with 
the  sugar-phosphate  backbone  of  the  DNA  strand  by  abstracting  an  H-atom  from  a  sugar 
molecule  to  form  water.  The  H  radical  can  also  abstract  H  molecules  from  the  sugars  on  the 
DNA  strand  to  form  H2.  Measured  reaction  rate  constants  of  the  OH  and  H  radicals  with 
the  sugar  moiety  of  the  DNA  strand  are  8  x  10^  M“^  s“*  and  7.8  x  10*  s“^  respectively 

[4]. 

The  reaction  of  OH  and  H  radicals  with  the  sugar-phosphate  DNA  backbone  is  one 
of  the  production  mechanisms  for  single  strand  breaks  (SSB),  which  are  scissions  of  one 
of  the  DNA  strands.  Scissions  of  both  strands  simultaneously  is  called  a  double  strand 
break  (DSB).  If  two  SSBs  occur  on  opposite  backbones  and  the  SSBs  lie  within  a  close 
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proximity  of  each  other,  the  two  SSBs  can  lead  to  a  DSB.  While  some  researchers  indicate 
that  the  maximum  separation  of  two  SSBs  that  can  become  a  DSB  is  six  to  ten  base  pairs 
[4],  others  claim  that  multiple  lesions  within  20  base  pairs  of  DNA  can  lead  to  a  DSB  [38]. 
Although  SSBs  are  often  easy  to  repair,  DSBs  are  difficult  to  repair  and  can  result  in  cell 
death,  mutations,  or  carcinogenisis. 

Indirect  action  has  also  been  shown  to  damage  the  structure  of  spore  layers.  Oxidizing 
agents,  to  include  H2O2,  appear  to  kill  spores  by  causing  damage  to  the  spore’s  inner 
membrane.  Identical  resistance  to  oxidizing  agents  was  seen  in  “spores  of  strains  with 
very  different  levels  of  unsaturated  fatty  acids”  [19].  Therefore  experimentalists  believe 
oxidative  damage  to  essential  proteins  in  the  inner  membrane  may  be  the  cause  of  spore 
death.  Researchers  suggest  that  spore  death  occurs  because  the  damaged  inner  membrane 
becomes  non-functional  or  ruptures  during  germination  [53].  Other  researchers  claim  that 
H2O2  reacts  with  molecules  in  the  spore  cortex  which  protects  the  core  from  H2O2  damage. 
However,  the  exact  method  in  which  this  would  inactivate  the  spore  is  unknown  [79]. 

Some  researchers  claim  lethal  doses  of  H2O2  inactivate  spores  by  damaging  its 
enzymes.  Experiments  have  shown  that  enzymes  necessary  for  germination  and  outgrowth 
have  experienced  lethal  damage  following  a  spore’s  treatment  with  H2O2  [73].  While  there 
is  no  consensus  on  the  precise  target  of  H2O2  reactions,  all  researchers  agree  that  H2O2 
does  not  inactivate  a  spore  by  damaging  its  DNA.  OH  has  been  shown  to  cause  damage 
to  eukaryotic  proteins  following  ionizing  irradiation.  And  while  this  damage  alone  was 
not  proven  to  be  lethal,  researchers  believe  that  ionizing  radiation-induced  protein  damage 
combined  with  DNA  breakage  may  result  in  the  inactivation  of  cellular  activities  necessary 
for  DNA  repair  [51].  Repair  of  DNA  damage  will  be  detailed  further  in  Chapter  4. 

3.6.2  UV  Radiation  Damage. 

The  primary  damage  mechanism  of  UV  radiation  is  direct  action  near  the  absorption 
site.  UV  radiation  produces  a  wide  range  of  DNA  damage  to  include  dimers  and  spore 
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photoproduct  (SP),  and  in  rare  cases,  SSBs  and  DSBs.  “Cyclobutane  pyrimidine  dimers 
constitute  the  major  DNA  photoproducts  upon  exposure  to  UV-B  light”  while  SP  is 
the  main  damage  mechanism  experienced  by  the  DNA  after  UV-C  irradiation  [61,  81]. 
However,  both  UV-B  and  UV-C  irradiation  of  bacteria  will  cause  dimers  and  SPs  to  form 
within  the  spore’s  DNA. 

DNA  bases,  primarily  the  pyrimidine  bases,  directly  absorb  incident  UV-B  photons. 
A  cyclobutane  pyrimidine  dimer  (CPD)  is  the  connection  of  two  adjacent  pyrimidine 
bases  on  the  same  DNA  backbone  strand.  CPDs  occur  when  absorbed  photon  breaks 
the  hydrogen  bonds  of  complementary  base  pairs  and  these  pyrimidine  bases  bond  with 
adjacent  pyrimidine  bases  on  the  same  DNA  strand.  The  three  types  of  CPDs  are  thymine- 
thymine,  thymine-cytosine,  and  cytosine-thymine.  Figure  3.2  demonstrates  the  formation 
of  a  thymine-thymine  CPD.  Dimers  will  create  a  kink  in  the  DNA  backbone  which  will 
prevent  the  DNA  from  being  transcribed  which  is  potentially  lethal  [105]. 


Figure  3.2:  Formation  of  Thymine  Cyclobutane  Dimer.  Figure  reprinted  with  permission 
from  [81].  Copyright  2001  Elsevier  Science  B.V. 


Spore  photoproduct  is  also  created  after  a  DNA  base  absorbs  an  incident  photon.  SP 
is  formed  between  two  adjacent  thymine  bases  on  the  same  DNA  strand  when  the  two 
bases  bond  to  the  same  methyl  group.  See  Figure  3.3  for  a  visual  representation  of  the  SP 
structure.  Approximately  30  percent  of  UV  irradiated  DNA  thymine  is  able  to  be  converted 
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into  SR  Similar  to  dimers,  SP  also  causes  helix-distorting  lesions  which  if  not  repaired,  can 
prevent  the  DNA  from  replicating  [90] . 


Figure  3.3:  Spore  Photoproduct.  Figure  reprinted  with  permission  from  [90].  Copyright 
1996,  American  Society  for  Microbiology. 


Small  acid-soluble  spore  proteins  (SASPs)  of  the  Q://?-type,  so  named  after  the  two 
major  proteins  found  in  Bacillus  subtilis  ,  provide  a  significant  contribution  to  the  spore’s 
UV  radiation  resistance.  The  a/jS-type  SASPs  are  double-stranded  DNA  binding  proteins 
that  each  cover  approximately  five  base  pairs.  There  are  enough  SASPs  in  each  spore  to 
saturate  the  spore  chromosome.  By  binding  to  the  outside  of  the  DNA  helix,  the  a/yS- 
type  SASPs  protect  the  backbone  from  enzyme  and  chemical  attacks.  Research  has  shown 
that  mutant  Bacillus  spores  lacking  the  cr//3-type  SASPs  are  much  more  UV  sensitive  than 
wild-type  (existing  in  nature)  spores.  It  is  believed  that  due  to  their  binding  capability,  a/yS- 
type  SASPs  prevent  the  formation  of  dimers  during  UV  irradiation,  however,  UV-B  still 
produces  enough  thymine-thymine  dimers  to  cause  spore  death  [90].  The  cr/jS-type  SASPs 
are  degraded  during  spore  germination,  and  if  not,  “their  continued  presence  interferes 
significantly  with  subsequent  development... by  blocking  DNA  transcription”  [2].  Thus,  the 
a/jS-type  SASPs  are  the  reason  that  UV  irradiated  dormant  spores  form  SPs  while  vegetative 
cells  mainly  experience  CPD  damage  following  UV  irradiation.  Now  that  radiation  damage 
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mechanisms  have  been  examined,  the  mathematical  modeling  of  radiation  damage  in  the 
literature  will  be  explored. 

3.7  Literature  Review 

Past  research  on  mathematical  modeling  of  ionizing  radiation  damage  has  focused 
on  the  effects  of  charged  particles  in  their  tracks.  The  charged  particles  create  radiation- 
induced  products  and  the  species  created  diffuse  and  react  with  the  matter  in  their  track. 
Thus  the  majority  of  ionizing  radiation  modeling  within  a  particle’s  track  has  utilized 
variations  of  the  reaction-diffusion  equation.  The  basic  reaction-diffusion  equation  of  a 
species  with  itself  in  the  absence  of  scavengers,  or  reactants  that  react  extremely  quickly 
with  other  like  species,  is 

—  =  DAR  -  kR^ 
ot 

where 

R  =  concentration  of  the  reactant  species 
t  =  time 

D  =  diffusion  coefficient 

k  =  rate  constant  for  a  reaction  of  the  species  with  itself. 

Many  one-radical  models  have  been  constructed  in  past  research  and  some  include 
the  existence  of  a  scavenger.  However,  models  that  only  consider  the  reaction-diffusion 
of  one  radical  are  chemically  unrealistic.  The  different  rate  and  diffusion  coefficients  are 
each  represented  by  one  value  and  many  models  did  not  recognize  the  difference  between 
radical  recombinations  and  other  chemical  reactions  which  misrepresents  chemical  yields 
[66]. 

Kuppermann  [66]  constructed  a  multiradical  diffusion  model  with  a  Gaussian  initial 
distribution  of  radicals.  The  model  contained  seven  reactive  species  and  three  scavengers. 
Kuppermann’s  model  was  solved  numerically  and  compared  to  experimental  data  with 
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varying  degrees  of  success.  Schwartz  [66]  also  constructed  and  numerically  solved  a 
multiradical  model  and  allowed  for  five  adjustable  parameters  related  to  initial  radical 
yields  and  spur  size.  The  adjustable  parameters  were  fit  to  sets  of  experimental  data. 

In  addition  to  the  number  of  radicals  and  scavengers  considered  in  a  model,  past 
ionizing  radiation  modeling  research  has  also  varied  regarding  the  geometry  of  the  track. 
Some  researchers  consider  a  cylindrical  track  while  others  claim  a  track  is  composed  of 
widely  spaced  spherical  spurs  which  can  be  thought  of  as  a  string  of  beads.  The  spherical 
spurs  strung  along  a  track  are  generated  by  low-LET  irradiation.  The  effect  of  high  LET,  as 
with  or-particles,  is  the  merging  of  spurs  to  form  a  cylindrical  track.  A  high  LET  cylindrical 
track  allows  for  nearly  complete  radical  recombination  while  a  spherical  spur  only  allows 
partial  recombination  [31]. 

Magee  and  Chatterjee  [31]  believe  that  the  chemical  product  yield  of  particles  should 
be  modeled  via  reaction-diffusion  equations  but  they  claim  that  the  track  should  be  modeled 
as  two  separate  entities.  Magee  and  Chatterjee  believe  the  track  consists  of  the  core  and  the 
penumbra  which  is  the  area  immediately  surrounding  the  core.  Within  the  core,  head-on 
collisions  occur  which  have  a  high  energy  loss  per  event  but  only  a  few  primary  events 
account  for  this  large  energy  loss.  Glancing  collisions  occur  within  the  penumbra  and  this 
type  of  collision  has  a  low  energy  loss  per  event  while  the  amount  of  events  is  relatively 
large. 

Stochastic  modeling  of  radiation  chemistry  has  included  Monte  Carlo  simulation, 
the  master  equation  approach,  and  the  independent  reaction  time  model.  Monte  Carlo 
simulations  model  the  reaction  and  diffusion  of  species  via  random  flight  simulations  and 
average  the  probabilities  of  reaction,  diffusion,  and  separation  over  a  large  number  of  trials. 
Eor  example,  if  a  particle  is  within  the  reaction  radius  of  another  particle,  it  is  assumed 
the  particles  will  react  with  a  given  probability.  Stochastic  approaches  which  allow  for 
fluctuation  in  the  number  of  reactant  species  have  been  suggested  as  a  valid  replacement  of 
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the  deterministic  differential  reaction-rate  equations.  Specifically  Gillespie  [36]  provided 
a  stochastic  simulation  algorithm  which  used  a  Monte  Carlo  procedure  to  model  collision 
probability  per  unit  time.  The  algorithm  considered  a  number  of  molecules  in  a  mixture  of 
chemical  species  which  can  interact  through  specified  chemical  reaction  channels  present 
at  an  initial  time.  The  model  predicted  what  the  molecular  population  levels  would  be  at 
any  later  time.  It  was  based  on  thermal  activation  energies  and  did  not  consider  diffusion. 
Monte  Carlo  simulation  is  computationally  intensive  and  has  not  yet  successfully  been 
compared  to  experimental  data  except  in  cases  with  very  few  reactants  [66] . 

The  master  equation  approach  “considers  the  state  of  a  spur  at  a  given  time  to  be 
composed  of  Ni  particles  of  species  f”  where  Ni  is  a  random  variable  [66].  This  method 
is  mathematically  similar  to  a  variation  of  the  diffusion  equation  and  it  can  be  easily 
generalized  for  multiradical  spurs.  However,  it  does  not  represent  the  various  reactions 
of  reactive  products  so  it  requires  a  relatively  simple  reaction  scheme.  The  independent 
reaction  time  model  treats  the  immediate  reaction  of  a  reactive  species  given  an  initial 
position  differently  than  its  intermediate  reactions  which  are  generated  by  the  inversion 
of  a  first  passage  time  problem.  The  model  assumes  that  pairwise  reaction  times  are 
independent  of  other  reactions.  Results  of  the  independent  reaction  time  model  have 
compared  positively  to  experimental  data  [66] . 

3.8  Mathematical  Damage  Model 

This  research  will  examine  radiation  damage  to  a  spore’s  DNA  and  enzymes  via  an 
illustrative  model  to  frame  all  of  the  important  indirect  radiation  damage  processes.  The 
model  described  has  relevance  to  sterilization  of  spores  or  eukaryotic  cells  adsorbed  on 
surfaces  in  an  ordinary  humid  environment.  First,  the  mathematical  model  of  damage  to  a 
spore’s  DNA  via  reaction-diffusion  partial  differential  equations  will  be  presented. 
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3.8.1  Reaction-Diffusion  ofROS  Within  the  Spore. 

A  demonstrative  scenario  will  be  explored  in  order  to  simulate  the  indirect  effects  of 
ionizing  radiation  in  a  spur  on  a  spore.  Consider  a  narrowly  collimated  source  of  radiation 
impinging  on  an  inert  solid  plate  having  a  thin  layer  of  adsorbed  water  of  thickness  h  and 
an  adsorbed  spore  (refer  to  Figure  3.4).  The  distance  between  the  radiation  source  and  the 
spore  is  considered  to  be  larger  than  the  spore  diameter  of  2r.  The  water  layer  is  thick 
enough  that  2r  <  h  and  yet  thin  enough  that  the  layer  can  be  assumed  to  be  uniformly 
saturated  with  dissolved  O2  adsorbed  from  the  air.  This  realizable  scenario  can  be  treated 
as  an  approximate  two  dimensional  model  of  radiation  induced  damage  to  a  spore  by  ROS. 


Adsorbed  Water 

Spor^ 

k 

f 

■ 

Inert  Solid 

Substrate 

Figure  3.4:  Illustration  of  Radiation,  Medium,  and  Spore  Scenario 


As  a  result  of  the  radiation  source,  several  reactive  species  are  generated  in  the 
neighborhood  of  the  radiation  path  through  the  O2  saturated  water.  These  reactive  species 
are  loosely  split  into  two  classes.  Class  1  will  contain  radicals  which  are  highly  reactive 
and  diffuse  rapidly  in  the  water  having  molar  concentration  Ri.  Ri  represents  the  class 
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that  reacts  quickly  with  the  O2  present  in  the  adsorbed  water  and  therefore  instantaneously 
transforms  into  other  species.  The  hydrated  electron,  e~g,  is  an  example  of  this  class  of 
reactants.  Class  2  will  contain  radicals  which  are  less  reactive  and  diffuse  slower  having 
concentration  R2.  R2  represents  the  class  that  does  not  react  as  fast  as  the  Ri  class  and 
therefore  does  not  get  transformed  as  quickly.  In  other  words,  if  k,  is  the  reaction  rate 
constant  for  species  Ri,  then  k2  ^  ky.  This  allows  the  class  of  R2  reactants  to  diffuse  to  the 
outer  boundary  of  the  spore.  HO2  and  examples  of  this  class  of  reactants. 

The  concentration  of  class  2  reactants,  R2,  can  be  determined  as  a  function  of  the 
dose  of  the  initial  stream  of  radiation.  Assume  this  concentration  is  uniformly  impressed 
on  the  outer  boundary  of  a  spherical  model  for  the  spore  as  it  would  be  if  the  spore  is 
located  far  enough  away  from  the  radiation  source  that  the  reactants  are  dispersed  uniformly 
throughout  the  substrate.  This  reactant  class  will  react  as  it  diffuses  into  the  spore.  The 
spore  is  assumed  to  have  uniform  material  properties  and  there  exist  no  permeability 
barriers  within  the  spore  that  would  interfere  with  the  diffusion  of  reactants  through  the 
spore.  Thus  the  spore  is  assumed  to  be  spherically  symmetric  (no  angular  dependence) 
and  the  diffusivity,  D3,  is  assumed  to  remain  constant  across  the  radius,  r,  of  the  spore. 
Finally,  it  is  assumed  that  initially  there  exist  no  radicals  within  the  spore  and  there  is  no 
flux  on  the  outer  boundary  of  the  spore.  The  process  will  be  modeled  as  a  partial  differential 
equation  with  boundary  and  initial  conditions  and  an  analytical  solution  will  be  presented. 
In  addition,  an  expression  for  the  average  reactant  concentration  within  the  spore’s  core 
will  be  discovered. 

The  partial  differential  equation  which  models  the  class  2  reactants  as  they  diffuse  into 
the  spore  is 

3R 

^  =  V  •  D2VR2  -  k3^2  (3.5) 

ot 

where 

Dj,  =  diffusivity  of  R2  into  the  spore 
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Riir,  t)  =  concentration  of  class  2  reactants 
r  =  spherical  distance 
t  =  time 

k-i  =  reaction  rate  constant. 


Note  that  R2  is  a  function  of  time,  t,  and  the  spherical  distance,  r.  The  term  k^Ri  represents 
the  rate  of  R2  reactions  within  the  spore.  In  other  words,  kj,  is  rate  at  which  class  2  reactants 
react  with  biological  compounds  within  the  spore  as  they  diffuse  throughout  the  spore.  The 
constant  diffusivity  and  lack  of  angular  dependence  allow  the  Laplacian  to  be  evaluated  in 
two  dimensions  and  Equation  3.5  becomes 


9Ri 


(3.6) 


Since  the  concentration  of  class  2  reactants  is  assumed  to  be  uniformly  impressed 
on  the  outer  boundary  of  the  spore,  the  initial  concentration  of  R2  outside  the  spore  is 
a  function  of  time  is  denoted  as  fit).  The  source  term  on  the  outer  boundary  gives  the 
following  boundary  condition 


Riiro,  f)  =  f{t)  where  represents  the  outer  spore  radius. 
Initially,  there  are  no  R2  radicals  in  the  spore,  so 


(3.7) 


^2(r,0)  =  0. 


(3.8) 


The  solution  must  be  finite  (particularly  at  r  =  0),  that  is, 

l^2(r,0l<«^.  (3.9) 

The  first  step  in  solving  Equation  3.6  is  applying  separation  of  variables.  Eetting 
R2{r,t)  =  X(r)T(t)  gives 


XT'  =  D3TX"  + 


ID^T  , 


-ksXT 


T 

Y 


D%  —  + 


r  X 


-k2  =  -n 


(3.10) 
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where  yu  is  a  separation  parameter.  The  solution  to  y  =  -yu  is 


T{t)  =  e 


(3.11) 


Upon  multiplication  of  Equation  3.10  by  it  becomes 


X"  +  -X'  +  - — —X  =  0 
r  D3 


which  written  in  self  adjoint  form  is 


\[r^X\r)\  +  =  0. 

H  D3 


(3.12) 


Let  X{r)  =  —  for  r  3^  0  so  that  r^X(r)  =  rJ(r).  In  addition, 


X'(r) 


rJ'ir)  -  J{r) 


Substituting  these  definitions  for  X{r)  and  X'{r)  into  Equation  3.12  produces 


Ip  pv  ^  jr  M/  ,  ^  “  ^3  J(r) 

—  [rJ  (r)  -  J(r)]  +  — - =  0 

Dt,  r 


which  reduces  to 


nr)  +  =  0. 


Now  let  A^in)  =  then  J"{r)  +  A^J(r)  =  0  and  J(r)  =  E  cos(/lr)  +  F  sin(/lr)  for  >  0. 


Di 


Eor  <  0,  only  the  trivial  solution  exists.  Therefore, 


03 


X{r) 


E  cos(Ar)  +  F  sin(/Ir) 


and 


X'(r)  = 


-EArsm{Ar)  +  FArco%{Ar)  -  E  cos(Ar)  -  F  sin(/lr) 


(3.13) 


It  has  been  assumed  that  the  spore  is  spherically  symmetric  therefore  at  the  center  of  the 


spore,  ^(0,0  =  0,  which  means  X'(0)T(t)  =  0.  Equivalently, 


lim  X'(r)  =  0. 

r^0+ 
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This  is  applied  to  Equation  3.13  and  leads  to  E  =  0  so  that 


X(r)  = 


F  sin(/lr) 


In  order  to  satisfy  the  homogeneous  boundary  condition,  X(ro)  =  0, 


x(.„) = = 0. 

rn 


(3.14) 


To  solve  this  equation,  note  that  F  =  0  will  produce  the  trivial  solution.  Therefore, 


sin(/lro)  =  0  so 


3  -  -10 

5  n  1,  Z, ... 


Thus  the  //„’s  are  the  eigenvalues  (recall  +  k^)  which  give  eigenfunctions,  X„’s, 


defined  as 


sin(d„r) 

X„(r)  =  Fn - ,  n  =  1, 2, ... 

r 


F„  is  found  by  normalizing  X„(r)  which  is  done  by  choosing  Fn  such  that 


JrTo  pin 

I  I  r^Xl  sin(0)d0d0dr 

0  Jo  Jo 


Jf'ro 

Fl 

0 

I 

0 


I  Sin^(/l„r)dr 


=  27rF„  [1 -cos(2d„r)]dr 


=  27rF„  To 


sin(2ro/l„) 


2nF^  ro-  —  sin(rod„)  cos(roT„) 


Equation  3.14  implies  sin(ro/I„)  =  0  which  means  that  1  =  InF^Vo  so  that 


Xnir) 


1 

Invo  r 
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Lemma  3.8.1.  The  eigenfunctions,  X„(r),  are  orthonormal  on  the  interval  (0,  r^)  with 
respect  to  the  inner  product  <  f,g  >=  4-n  r^f(r)g(r)dr  ( i.e.  <  >=  6mn)- 

Proof: 


<  Xn,X„r  >  =  47T  f  r^X„(r)X„(r)dr 
Jo 


,  I  /  -  ■  I  \  \  ^ 

47r  I  r- - r  sin  — r  sin  — r  dr 


r  2^ 

Jo  ^  2;rr, 

2  r'^"  .  Inn  \  .  ( mn  \  , 
=  —  I  sin  — r  sin  — r  dr 

ro  Jo  \ro  j  \ro  j 


if  n  =  m  then 


Xfi ,  > 


2  n  .  Jnn 
'o  Jo  \ 

7/ 


1  I 

1  -  cos  - r 


dr 


r„  - 


r^  sin  {Inn) 
Inn 


and  because  sin(2n;r)  =  0,  this  case  gives  <  X„,Xm  >=  1.  Ifn  4^  m  then 

2  r'^°  .  Inn  \  .  I  mn  \  , 

<Xn,Xyn>=—  I  Sin  — r  Sin  — r  dr  =  0 

fo  Jo  \^o  j  \  j 

therefore  combining  the  two  cases  gives 

<  Xfi ,  Xffi  >  —  5fYin  •  B 

Since  +  k^.  Equation  3.11  becomes 

Tn{t)  =  A„e“^7  n  =  1,2,  ... 

where  A„  is  an  arbitrary  constant.  Therefore,  the  solution  to  the  homogeneous  partial 
differential  equation  with  homogeneous  boundary  conditions  is  of  the  form 

^  ^  -fInFo  r 

n=\  n=\  ’  ^ 
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Next,  the  forcing  function  at  the  boundary,  7?2(ro,  t)  =  f{t),  should  be  considered.  Let 


R2{r,t)  =  f{t)  +  W{r,t) 


where  W{r,  t)  =  2^=1  ^n{r)Tn{t).  Note  that  w(ro,  t)  =  0  and 

/ST?  ^ 

=  fit,  d)  +  W'ir,  t)  =  fit)  +  ^  Xniffit). 

n=\ 

Recall  that  the  original  partial  differential  equation  (Equation  3.5)  was 

dR2 

-  V  •  D3VR2  -  ^3^2 
ot 


which  can  combine  with  Equation  3.16  to  give 


00  00  00 

fit)  +  ^  Xniffit)  =  -hfit)  -hYj  ^nifTnit)  +  ^  Tft)^  [r^Kif]  • 

n=\  n=\  n=l 

Using  Equation  3.12  for  the  last  term  in  Equation  3.17  gives 

00  00  cx> 

fit)  +  X„ir)fit)  =  -hfit)  -hYj  ^nifTnit)  +  fmlXfr) 

n=l  n=l  n=\ 


J]  [fit)  +  hfit)  +  D^XlTnit)]  Xnir)  =  -hfit)  -  fit). 

n=\ 

After  applying  Eemma  3.8.1, 

fit)  +  tlnTnit)  =  [-^3/(0  -  fit)]  <hXn> 
where,  using  integration  by  parts. 


4  rl 

= - —=  —  cosinn) 

^3/2 

_  ^ '  o 

n 
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(3.15) 


(3.16) 


(3.17) 


(3.18) 


Multiply  both  sides  of  Equation  3.18  by  the  integrating  factor  e^"'  produces 


e^"'[r;(o  +yu„r„(o]  =  -  [e^"^r„(o] 

=  en-W)-/'(0]<l,^„> 

Integrating  both  sides  of  the  equation  with  respect  to  t  gives 

e^-'r^CO  -  r„(0)  =<  1,  >  r  -  e^"7'(T)]  dr 

Jo 

Integration  by  parts  applied  to  the  right  hand  side  of  the  equation  produces 

-  r„(0)  =<  1,X„  >  -e^'-%f(T)dT  -  [e^"V(r)]o  +  ^  yu„e^ 

=<  1,X„  >  I  7'e""7(r)dr  -  e""7(0  +  /(0)|  • 

Initially  (at  time  t  =  0)  there  are  no  radicals  in  the  spore,  i.e.  R2{r,  0)  = 
Equation  3.15  becomes 

oo 

R2(r,  0)  =  0  =/(0)  +  X„(r)TM 

n=i 

00 

-m  =J]Xn(r)Tn(0). 

n=\ 

Testing  with  and  applying  Eemma  3.8.1  results  in 


r„(0)  =  -  <  i,x„  >  /(O) 

Using  this  initial  condition,  Equation  3.19  becomes 

e^"7„(0  =<  1,X„  >  d^e^"7(T)dT-e^"7(oj 


or 


Tnit)  =  -<l,Xn>  f(t)+  <l,Xn>  f  T^e^"("-7(T)dT. 

Jo 


’7(T)dT| 

(3.19) 
0.  Therefore, 
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Now 


R2{r,t)  =  f{t)  +  Y^Xn{r)Tn{t) 

n=\ 


=  fit)  +  <  1,X„  >  f(t)+  <  l,X„  > 

n=l  ^ 

OO 

=  yiXW<l,X„>  e-^'-(^-y(T)dT 

^  -'0 


/ 


^y^(r-0f(T)dT 


(3.20) 


as 

00 

J^X„(r)<l,X„>=l. 

n=l 

The  spore’s  DNA  is  located  within  its  core.  Therefore,  to  determine  the  amount  of 
reactants  that  can  damage  the  spore’s  DNA,  the  average  concentration  of  class  2  reactants 
within  the  core  must  be  calculated.  Let  denote  the  average  concentration  of  class 
2  reactant  within  the  core  of  the  spore.  In  order  to  calculate  the  average  concentration, 
the  total  R2  concentration  in  the  core  is  found  first  by  integrating  with  respect  to  spherical 
coordinates  over  the  core  volume.  Then  the  total  concentration  is  divided  by  the  core 
volume  (Vcore  =  where  rc  denotes  the  radius  of  the  spore’s  core. 


1 

^2(0  =  i)r^  sin(0)d0d0dr 

^core  ^0  ^0  ^0 


Jf>Yc 

R2(r,t)r^  I  I  sin(^)d0d0dr 
0  Jo  Jo 


47rr^ 


3  r'"^ 

=  —  I  R2(r,  t)r^dr. 

Jo 

Substituting  Equation  3.20  into  Equation  3.21  gives 

R^^it)  =  e-'^"('-y(T)dr|r2dr 

=  ^  E  {X  {X  ^  >  ^"drj 

=  —  y  [-2cos(n;r)]  I  e~^"‘^^“^V(T)dT  I  /l„r  sin(/l„r)dr 

A  ^  Jo  Jo 

=  ^  ^[-2cos(n;r)]  sin(T„rc)  -  rcCOs(T„rc)|  e“^"('“y(r)dr 


(3.21) 


(3.22) 
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which  reduces  to 


Rl{t)  =  yKn  e-^"('-^V(r)dr  (3.23) 

^  -'o 

where 

Kn  =  \Al<l,Xn>  \\n(rVAr  (3.24) 

Jo 

=  (-1)”+'  ^  sin(d„rc)  -  r,  cos(d„rc)| .  (3.25) 

The  parameters  utilized  within  the  reaction-diffusion  model  are  included  in  Table 
3.5.  The  initial  concentration  of  R2  on  the  outer  boundary  of  the  spore,  denoted  by  f(t), 
is  developed  from  experimental  values  and  models  of  species’  reaction-diffusion  within 
an  ionizing  particle’s  spur.  The  collimated  source  is  assumed  to  be  low-LET  radiation 
applied  via  a  pulse  (1  nanosecond  (ns)  pulse  of  20  Mev  electrons).  Between  10“'^  and 
10“^^  seconds  after  irradiation  of  an  aqueous  solution,  ions  and  excited  states  are  formed. 
Chemical  reactions  begin  to  take  place  around  10“^^  seconds  and  radicals  within  a  spur 
reach  their  initial  yield,  or  highest  concentration,  at  around  10“  seconds  as  shown  in 
Figure  3.5  [44,  66].  At  this  time,  the  chemical  yields  for  Off,  ff,  and  ff2  are  5.9,  0.62,  and 
0.15  molecules  per  100  eV  respectively  [31].  Off  decays  at  a  fast  rate  after  it  reaches  its 
initial  yield  (i.e.  the  Off  yield  decreases  by  30  percent  within  3  ns  where  1  ns  is  equivalent 
to  1  X  10“^  seconds)  which  is  also  depicted  in  Figure  3.5  [31,  44].  Note  that  the  time  scales 
for  the  horizontal  axis  in  Figure  3.5  are  different. 

Table  3.5:  Parameters  of  the  Reaction-Diffusion  Model 


Parameter 

Value 

Source 

ro 

.58  (//m) 

[106] 

rc 

.4  (jum) 

[106] 

h 

.5  (nM-i  s-^) 

[4,  30] 

D3 

2  X  10^  (yum^ 

[31] 
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Time  (s)  Time  (s)  x'<>® 

(a)  Production  of  Class  2  Radicals  (b)  Decay  of  Class  2  Radicals 

Figure  3.5:  Concentration  of  R2  on  outer  boundary  of  spore 


The  parameters  in  Table  3.5  and  the  initial  concentration  of  R2  on  the  outer  boundary 
of  the  spore  shown  in  Figure  3.5  are  utilized  to  find  the  average  concentration  of  class  2 
reactants  in  the  spore  core  given  by  Equation  3.23.  Figure  3.6  demonstrates  the  evolution 
of  Rj  over  time.  Note  that  the  core  concentration  of  reactants  reaches  its  maximum  value  at 
approximately  10  ns  after  the  application  of  the  low-LET  radiation  pulse  and  then  decays 
rapidly  due  to  the  decrease  of  R2  reactants  on  the  spore  boundary.  Radicals  in  the  core  will 
damage  the  spore’s  DNA  and  enzymes  which  will  be  discussed  in  the  next  two  sections. 

3.8.2  DNA  Damage. 

The  rate  of  DNA  damage  caused  by  ROS  is  modeled  by 

-  ^  =  kdR\{t){D]{t)  (3.26) 

at 

where 

[D](f)  =  concentration  of  DNA  in  the  spore  (nM) 
t  =  time  (s) 

kd  =  rate  coefficient  of  ROS  interaction  with  DNA  ^nM“^s“*j 


42 


0  0.1  02  0.3  0.4  0.5  0.6  0.7  0.8  0.9  1 


Time  (s)  xio-^ 

Figure  3.6:  Average  Coneentration  of  R2  Reaetants  in  Spore  Core 


and  Ry<t)  is  the  average  eoneentration  of  class  2  reactants  in  the  core  defined  by  Equation 
3.23.  Equation  3.26  characterizes  the  rate  of  change  in  the  DNA  content  of  a  spore.  By 
integrating  with  respect  to  time,  the  solution  to  the  rate  equation  is 

[D]a)  =  {Dhfoit)  (3.27) 

with  the  fraction  of  DNA  remaining  after  some  radiation  exposure  time,  t,  defined  as 
fait)  =  where 

R(t)  =  f  R^(T)dT.  (3.28) 

Jo 

At  time  t  =  0,  the  initial  concentration  of  DNA,  [D]t,  is  determined  by  a  spore’s  total 
DNA  concentration  prior  to  radiolysis  which  is  given  by  the  value  2.87  x  10^  nM  [11].  The 
initial  DNA  concentration  in  the  core  is  based  on  the  concentration  of  DNA  in  the  E.  coli 
cell.  It  is  utilized  for  the  purpose  of  uniformity  because  while  the  concentration  of  DNA 
in  the  B.a.  spore  is  known  (3.23  x  10’  nM  [8]),  repair  rate  parameters  used  in  Chapter  4 
are  more  readily  available  for  the  E.  coli  cell.  The  parameter,  kd,  is  given  by  the  reaction 
rate  of  the  hydroxyl  radical  with  the  sugar  moiety  of  the  DNA  strand  (i.e.  0.8  nM“^  s“  J 
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[4].  Figure  3.7  illustrates  the  fraetion  of  undamaged  DNA,  /d,  remaining  in  the  spore  over 
time.  As  the  radicals  in  the  core  react  with  and  damage  the  spore’s  DNA,  the  percentage 
of  viable  DNA  decreases.  The  irradiation  via  a  low-LET  pulse  (1  ns  pulse  of  20  Mev 
electrons)  damages  approximately  10  percent  (exactly  9.63  percent)  of  the  DNA  according 
to  the  model  presented  here. 


Figure  3.7:  Fraction  of  Undamaged  DNA  in  Spore  After  Irradiation 


3.8.3  Enzyme  Damage. 

In  addition  to  DNA  damage,  damage  is  inflicted  on  the  spore’s  proteins  during 
radiation  exposure.  As  discussed  in  Section  3.6.1,  ROS  have  been  proven  to  damage  the 
spore’s  enzymes  which  are  proteins  that  catalyze  chemical  reactions.  DNA  repair  proteins 
existing  within  the  spore  are  necessary  to  repair  damage  accumulated  within  the  dormant 
spore  and  thus  are  extremely  important  to  spore  outgrowth  [42].  Since  “accumulation  of 
oxidative  damage  to  proteins  diminishes  their  catalytic  activities  and  interactions”  [22], 
some  sources  claim  the  death  of  a  cell  is  more  highly  correlated  to  protein  damage  than  to 
DNA  damage  [50].  Others  go  as  far  as  to  state  that  the  spore’s  proteins  are  the  “principal 
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target  of  the  biological  action”  of  ionizing  radiation  [23].  Therefore,  DNA  must  be  repaired 
before  a  viable  cell  is  produced  during  germination  and  there  must  be  a  suitable  amount  of 
repair  enzymes  that  survive  the  radiation  exposure  to  accomplish  this.  Since  ROS  reacting 
with  DNA  and  enzymes  could  decimate  these  molecules  to  such  an  extent  that  they  can  no 
longer  be  repaired,  an  enzyme  damage  model  is  presented. 

There  are  six  enzymes  involved  in  nucleotide  excision  repair  (NER)  which  is  the 
primary  mechanism  of  radiation  repair.  NER  and  the  role  each  enzyme  plays  in  the  repair 
process  will  be  described  in  Section  4.3.  Eet  e{t)  be  the  vector  of  enzyme  concentration  for 
a  given  spore  within  the  population  at  time  t  where  e{t)  =  [exit), . . .  ■  The  rate  of 

enzyme  damage  caused  by  reactions  with  ROS  is  represented  by 

-%=  K,R^{t)e{t)  (3.29) 

at 

where  R^it)  is  defined  by  Equation  3.23  and  Ke  is  a  diagonal  matrix  of  rate  coefficients  for 
reactant  interaction  with  each  enzyme  (nM“^s“^j  i.e.  Ke=d.mg{ke^, . . .  ,ke^].  This  equation 
characterizes  the  rate  of  change  in  viable  enzyme  activity  of  a  spore.  By  integrating  with 
respect  to  time.  Equation  3.29  becomes 

e{t)  =  meo  (3.30) 


where 


m  = 


is  the  diagonal  matrix  for  the  fraction  of  remaining  enzymes  after  irradiation.  The  vector 
^0  =  [^io>  •  •  •  >  ^6o]^  is  the  initial  concentration  and  R(t)  is  defined  in  Equation  3.28.  Eor  this 
work,  the  rate  coefficients,  k^,  are  assumed  to  have  the  same  value  for  each  enzyme  and 
are  further  assumed  to  be  identical  to  kd  which  is  the  reaction  rate  of  the  OH  radical  with 
DNA. 

At  time  t  =  0,  prior  to  radiation  exposure,  the  initial  enzyme  concentration,  cq 
is  determined  by  the  spore  chosen  from  a  population  with  a  certain  enzyme  content 
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density  function  based  on  various  preparation  parameters.  “Spores  prepared  at  different 
temperatures  have  a  number  of  significant  differences  in  properties”  [35].  In  addition, 
spore  preparation  temperature  effects  the  core  water  content  which  in  turn  influences  a 
spore’s  resistance  to  radiation  damage  [59].  Therefore,  this  research  will  consider  that 
a  population  of  radiation  treated  spores  has  a  certain  initial  ‘fitness’  probability  density 
function  (PDF)  based  on  viability  of  enzyme  capacity  to  repair  DNA,  and  by  choosing  one 
spore  from  the  sample,  we  are  provided  a  specific  enzyme  concentration.  At  time  t  =  0, 
Cq,  is  randomly  drawn  from  this  continuous  PDF  which  has  a  mean  of  fieo  and  a  variance 
denoted  by  Note  that  the  initial  enzyme  ‘fitness’  is  assumed  to  be  independent  of  the 
fraction  of  remaining  enzymes  after  irradiation. 

The  mean  and  variance  of  a  continuous  random  variable  can  be  defined  in  terms  of 
expectations.  Thus 

lieo  =  E[eo]  (3.31) 

and 

o-l  =  V[eo]=E[{eo-fie,f].  (3.32) 

Observe  that  within  the  variance  expression,  the  square  of  eo  -  Heo  is  done  component 
wise.  Expectation  is  a  linear  operator  so  given  constants  a  and  {3,  E[QrX+;S]=aE[X]+jS  for 
a  random  variable  X  [17].  Since  eo  is  a  random  variable,  e(t)  is  also  a  random  variable 
drawn  from  a  PDE  with  mean,  /ig,  and  variance,  cr^.  Using  the  linear  operator  properties 
and  Equation  3.30,  the  mean  of  e(t)  is  defined  as 

fig  =  E[e(t)] 

=  E  [fg(t)eo] 

=  fe(t)E[eo] 

~  fe(3)flgQ. 
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In  addition,  the  variance  of  e{t)  is  given  by 


(rl  =  E  [(e(0  - 

=  {Ut)fE[{eo-Eeof]  (3.33) 


Note  that  the  development  of  Equation  3.33  is  possible  because  fe{t)  is  a  diagonal  matrix  of 
constant  values.  So  the  variance  of  the  e{t)  PDF,  cr^,  is  scaled  by  a  matrix  factor  of 

The  six  NER  enzymes  are  UvrA,  UvrB,  UvrC,  helicase  II,  polymerase  I,  and  ligase. 
Values  found  in  the  literature  for  the  number  of  UvrB  molecules  in  an  E.  coli  cell  range 
from  approximately  200  to  250  and  approximately  20  to  25  copies  of  the  UvrA  enzyme  are 
reported  to  be  in  each  cell  [100,  103].  Also,  researchers  have  shown  that  approximately 
2400  to  3000  molecules  of  the  helicase  II  enzyme  exist  in  each  E.  coli  cell  [7,  76].  Finally, 
479  copies  of  polymerase  I  and  226  copies  of  ligase  in  the  E.  coli  cell  were  counted 
during  an  experiment  while  other  literature  reports  there  are  approximately  400  and  200 
molecules  per  E.  coli  cell  of  polymerase  I  and  ligase,  respectively  [101].  Therefore, 
the  spore  population’s  initial  concentration  for  each  of  the  six  enzymes  is  modeled  via 
the  Uniform  PDF  with  a  range  of  20  percent  about  their  mean  values  (i.e.  the  range  is 
plus/minus  10  percent  from  the  mean).  The  state  of  each  of  the  six  enzymes  is  assumed  to 
be  independently  distributed.  The  value  used  for  the  mean  of  each  enzyme  of  the  E.  coli 
cell  is  outlined  in  Table  3.6.  Note  that  both  E.  coli  and  the  Bacillus  genus  utilize  the  same 
six  NER  enzymes  however  the  number  of  copies  for  all  of  the  enzymes  in  the  E.  coli  cell  is 
reported  in  the  literature  and  for  consistency,  is  used  here. 

The  Uniform  PDF  is  represented  by 


(pxix)  =  I 


1 

b-a^ 

0, 


a  <  X  <b 
otherwise 


(3.34) 
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Table  3.6:  Enzyme  Data  for  E.  coli 


Enzyme 

Copies/Cell 

Source 

Concentration 

A 

23 

[100,  103] 

38.2  nM 

B 

225 

[100,  103] 

373.6  nM 

C 

10 

[100,  103] 

16.6  nM 

H 

2750 

[7,  76] 

4566.5  nM 

P 

440 

[101] 

730.6  nM 

F 

200 

[101] 

332.1  nM 

where  b  -  a  represents  the  range  of  X,  the  continuous  random  variable.  The  Uniform  PDF 
has  mean,  £'[X]  =  ^  and  variance,  y[X]  =  [17].  Figure  3.8  is  an  illustration  of  the 

density  function  for  the  population’s  UvrB  content  prior  to  and  after  radiation  exposure. 
The  red  line  represents  the  PDF  of  the  population’s  undamaged  UvrB  concentration  which 
has  a  mean  of  373.6  nM  and  a  range  of  336.2  to  411.0  nM.  The  black  PDF  depicts  the 
distribution  of  the  UvrB  enzyme  when  it  has  been  degraded  to  90  percent  of  its  original  state 
(i.e.  its  mean  is  now  336.2  nM).  The  blue  line  illustrates  the  UvrB  enzyme  at  70  percent  of 
its  initial  concentration  (i.e.  the  PDF  mean  is  261.5  nM).  Notice  that  the  variance  of  each 
PDF  decreases  as  the  damage  inflicted  to  the  enzyme  increases  so  that  as  the  activity  level 
decreases,  the  population  begins  to  behave  more  like  the  population  mean. 

Assuming  the  six  NER  enzymes  function  independently,  their  joint  density  function 
of  enzyme  concentration  values  is  given  by 

(  6 


6 

(p  {e{t))  =  Y\  (pi{ei{t))  =  < 

i=\ 


n  bi-at  ’ 


i=\ 

0. 


where  e  represents  the  continuous  random  variables. 


at  <  ei  <  bi 


(3.35) 


otherwise. 

Prior  to  any  radiation  damage  (i.e. 


48 


0.02 


0.018 


0.016 


0.014 


^  0.012 
S 

ro  0.01 
P 

>  0.008 


0.006 


0.004 


0.002 


0 

200  250  300  350  400  450  500 

UvrB  Concentration  (nM) 

Figure  3.8:  PDFs  for  Undamaged  and  Damaged  UvrB  Enzyme 


'  Undamaged  PDF 
'  10%  Damage 
'  30%  Damage 


when  t  =  0),  the  Uniform  PDFs  of  the  enzyme  eoneentrations  have  a  20  pereent  span  about 
their  mean,  ■  Thus 


ai  =  fe(t)  (0.9yUg,  J, 

bi  = 

and  the  range  of  the  eontinuous  random  variable,  e,  ,  is 


bi  -  at  =  0.2 


Therefore  the  joint  density  funetion  of  enzyme  eoneentration  values  given  in  Equation  3.35 


beeomes 


(p  ie(t))  =  < 


1 

1=1  ^ 


0, 


0.9fe{t)n,.^  <  et  < 
Otherwise. 


(3.36) 
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The  domain  of  the  joint  density  function  represents  a  ‘hypercube’  of  possible  enzyme 
combinations  that  are  available  to  repair  damaged  DNA.  This  ‘hypercube’  will  be  explored 
further  when  building  the  population’s  probability  of  survival  model  in  Chapter  5. 

3.9  Summary 

This  chapter  examined  the  indirect  damage  to  a  spore  due  to  ionizing  radiation.  The 
primary  insult  mechanism  considered  was  indirect  damage  from  the  production  of  reactive 
oxygen  species  (ROS).  These  reactants  diffuse  into  the  spore  core  and  not  only  react  with 
the  DNA  but  also  with  the  enzymes  necessary  for  repair  of  the  DNA.  The  evolution  of 
the  reactive  oxygen  species’  within  the  core  was  determined  based  on  a  reaction-diffusion 
equation.  Using  the  average  core  concentration  of  ROS,  a  damage  model  for  both  DNA 
and  essential  enzymes  was  considered.  However,  DNA  repair  mechanisms  are  extremely 
efficient  thus  the  spore’s  repair  of  radiation  damage  must  be  taken  into  account  before 
considering  its  probability  of  survival.  DNA  repair  methods  and  a  mathematical  repair 
model  will  be  outlined  in  the  next  chapter. 
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IV.  Repair  Model 


A  spore’s  damaged  DNA  is  repaired  via  an  intrieate  system  involving  interaetion 
of  multiple  moleeular  proeesses.  There  are  five  outcomes  from  DNA  repair,  the  first  of 
which  is  that  the  damage  is  repaired  correctly.  The  second  occurs  when  DNA  damage  is 
repaired  incorrectly  but  the  damage  did  not  take  place  in  a  critical  location  necessary  for 
replication  or  transcription,  so  the  spore  may  experience  normal  outgrowth.  For  example 
a  mutation,  which  is  a  change  in  the  DNA’s  base  sequence,  may  occur  to  the  third  base 
which  encodes  a  codon.  However,  due  to  redundancy  in  the  genetic  code,  this  may  not 
effect  the  identification  of  its  corresponding  amino  acid.  The  third  result  of  DNA  repair  is 
that  the  damage  is  repaired  incorrectly  and  it  does  lie  in  a  critical  location.  In  this  case,  the 
DNA  may  replicate  a  mutation  which  will  then  be  passed  onto  future  generations  and  only 
a  fraction  of  cells  survive  this  mutagenic  process.  Or  the  incorrectly  repaired  damage  may 
prevent  the  spore  from  replicating  and  subsequent  outgrowth  which  will  lead  to  spore  death. 
Finally,  the  DNA  may  not  be  repaired  at  all  because  sufficient  proteins  were  not  available 
within  the  spore  for  the  repair  or  replication  process.  This  scenario  will  be  discussed  further 
in  Section  4.5. 

As  discussed  in  the  previous  chapter,  damage  to  a  spore’s  DNA  includes,  but  is  not 
limited  to,  base  modifications,  SSBs,  and  DSBs.  As  the  range  of  radiation  damage  to 
DNA  is  very  broad,  so  too  are  the  DNA  repair  mechanisms.  The  four  main  DNA  repair 
pathways,  all  of  which  incorporate  a  series  of  enzymatic  processes,  are  base  excision  repair 
(BER),  mismatch  repair,  nucleotide  excision  repair  (NER),  and  DSB  repair.  This  chapter 
will  describe  these  four  repair  pathways  as  well  mathematically  model  the  NER  process 
because  it  repairs  the  primary  DNA  damage  mechanisms  of  ionizing  radiation. 
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4.1  Base  Excision  Repair  (BER) 

The  BER  mechanism  corrects  base  modifications  such  as  depurination  and  deamina¬ 
tion  in  which  the  DNA  base  is  chemically  altered,  destroying  its  identity.  The  BER  process 
is  initiated  by  proteins  which  recognize  a  damaged  base  such  as  DNA  glycosylase.  Then 
the  damaged  base  is  removed  and  the  DNA  backbone  is  nicked  on  both  sides  of  the  cleaved 
base.  These  two  steps  are  completed  by  one  or  more  of  the  lyase  enzymes.  Next,  the  correct 
complementary  base  pair  is  filled  in  by  DNA  polymerase  and  finally,  the  DNA  backbone  is 
resealed  via  the  ligase  enzyme  [39,  105]. 

4.2  Mismatch  Repair 

Mismatch  repair  corrects  errors  caused  during  DNA  replication.  Because  mismatch 
repair  occurs  after  replication,  the  damage  detection  enzymes  must  be  able  to  identify 
which  DNA  strand  contains  the  replication  damage  and  which  one  is  the  original  “correct” 
strand.  It  is  possible  for  mismatch  repair  to  misidentify  the  two  strands  and  thus  “fix”  the 
wrong  strand.  This  causes  a  mutation,  or  change,  in  the  DNA  which  may  be  replicated  in 
subsequent  daughter  chromosomes  [105].  With  the  exception  of  occurring  after  replication, 
the  mismatch  repair  process  is  very  similar  to  NER. 

4.3  Nucleotide  Excision  Repair  (NER) 

The  NER  process  repairs  a  section  of  damaged  DNA  and  it  is  responsible  for  the 
removal  of  a  broad  group  of  structurally  unrelated  DNA  lesions  to  include  SSBs  and  the 
interactions  of  damaged  bases.  Disruption  of  DNA  base  pairing,  DNA  bendability,  and 
the  presence  of  chemical  modifications  within  the  DNA  are  all  fixed  via  NER  [37].  In 
addition,  bacterial  spores  use  NER  to  repair  dimers  [105]  and  spore  photoproduct  (SP)  and 
the  repair  process  is  “efficient  and  relatively  error  free”  [68].  Recall  from  Section  3.6.2 
that  dimers  and  SP  are  the  most  likely  types  of  damage  produced  by  of  UV  radiation.  The 
primary  DNA  damage  mechanisms  due  to  ionizing  radiation  are  base  modifications  and 
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SSBs  so  this  research  will  mathematically  model  the  NER  mechanism.  As  such,  NER  will 
be  discussed  in  detail. 

The  NER  process  uses  two  broad  stages  to  repair  damage.  The  first  involves 
recognition  and  incision  of  the  location  with  the  DNA  strand  containing  the  lesion  and 
the  second  consists  of  restoring  the  original,  undamaged  DNA  [6].  As  the  steps  of  NER  are 
described,  refer  to  Eigure  4. 1  for  their  pictorial  representation. 

The  first  step  in  the  NER  process  is  damage  recognition  which,  in  bacteria,  is  done 
by  the  proteins  UvrA  and  UvrB.  The  UvrA  protein  is  responsible  for  probing  the  DNA 
for  damage  and  then  signaling  the  location  of  possible  damage  to  the  UvrB  protein.  The 
presence  of  a  lesion  is  then  verified  by  the  UvrB  protein.  The  UvrA  and  UvrB  proteins 
form  a  UvrA2B2  or  a  UvrA2Bi  complex  (see  ‘A2B1-DNA  complex’  in  Eigure  4.1).  The 
UvrA2B  complex  is  a  member  of  the  helicase  super-family  and  therefore  also  unwinds  and 
separates  the  DNA  at  the  damage  site  and  signals  the  location  of  the  damage  to  the  UvrC 
protein.  UvrC  nicks  the  DNA  backbone  on  both  sides  of  the  damage  (see  ‘5’  incision  by 
UvrC’  in  Eigure  4.1).  The  UvrA2B  complex  acting  in  coordination  with  the  UvrC  protein 
is  often  referred  to  as  the  UvrABC  endonuclease  [37,  43,  75]. 

The  next  step  of  NER  is  completed  by  the  protein  UvrD  (also  called  helicase  II)  which 
displaces  the  damaged  DNA  strand.  Then  the  enzyme,  polymerase  I  (or  Pol  I)  replaces  the 
excised  DNA  by  laying  down  the  correct  complementary  base  pairs.  Einally,  DNA  ligase 
seals  the  nicks  of  the  DNA  backbone  strand  by  connecting  the  newly  corrected  section  of 
the  backbone  to  the  old  backbone  (see  ‘Repair  synthesis  and  ligation’  in  Eigure  4.1)  [105]. 

4.4  Double  Strand  Break  Repair  (DSB) 

DSBs  are  repaired  by  one  of  two  methods;  recombination  repair  or  non-homologous 
end-joining.  Bacteria  can  only  initiate  recombination  repair  after  DNA  replication  because 
it  requires  two  copies  of  the  DNA.  Recombination  repair  involves  copying  DNA  sections 
of  the  daughter  chromosome  and  using  them  to  replace  the  sections  of  DNA  containing  the 
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Repair  synthesis  - •  • 

and  ligation  - 


Figure  4.1:  Nucleotide  Excision  Repair.  Figure  reprinted  with  permission.  This  image  was 
published  in  [75].  Copyright  Elsevier  (1999). 


DSB.  While  recombination  repair  is  the  preferred  method  of  repair  for  DSBs,  there  exists 
significant  complexity  in  identifying  correct  corresponding  sequences  and  facilitating  their 
cross  over  to  the  damaged  strands.  Non-homologous  end-joining  (NHEJ)  directly  rejoins 
the  DNA  chains.  This  process  aligns  the  broken  ends  of  the  DSB,  trims  any  frayed  ends, 
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and  ligates  the  strands.  Recombination  repair  is  the  preferred  mechanism  for  repairing 
DSBs  because  NHEJ  is  very  error  prone  [105]. 

4.5  Literature  Review 

To  date,  mathematical  modeling  of  DNA  repair  has  almost  solely  focused  on  human 
cells.  In  addition,  much  of  the  research  has  concentrated  on  modeling  the  DNA  repair 
processes  as  a  whole  by  grouping  together  all  the  DNA  repair  mechanisms  into  one 
model  and  conducting  some  sort  of  data  fitting.  Foray  et  al.  [33]  proposed  a  stochastic 
and  multicomponent  model  for  repair  of  human  cells  based  on  “the  existence  of  two 
subpopulations  of  DNA  damage  with  a  fast  and  a  slow  constant  repair  half-time”  [34]. 
They  showed  that  the  Gamma  distribution  most  accurately  fit  this  hypotheses  as  compared 
to  laboratory  measurements.  Gastaldo  et  al.  [34]  upheld  this  DNA  repair  rate  interpretation 
and  applied  it  to  a  combined  DNA  damage  induction  and  repair  model.  Theoretical  and 
experimental  data  was  again  analyzed  however  this  theory  also  did  not  distinguish  between 
the  various  repair  mechanisms  required  to  fix  the  multiple  types  of  DNA  damage  such  as 
base  modifications,  SSBs,  and  DSBs  incurred  after  ionizing  irradiation. 

The  first  mathematical  model  that  closely  examined  a  single  repair  mechanism 
was  proposed  by  Sokhansanj  et  al.  [94].  They  mathematically  modeled  the  human 
BER  pathway  as  a  multi-step,  sequential  process  using  Michaelis-Menton  kinetics  and 
quantitative  data  on  single  proteins.  In  addition,  the  sensitivity  of  the  pathway  to  altered 
enzyme  kinetics  was  analyzed.  Crooke  et  al.  [21]  also  presented  a  model  of  the  human 
BER  process  and  incorporated  the  damage  stage  with  the  repair  stage  in  the  same  model. 
Their  model  included  a  set  of  differential  equations  representing  the  sequence  of  enzymatic 
reactions  in  both  the  damage  and  repair  pathways  using  Michaelis-Menton  kinetics.  This 
model  assumed  a  certain  initial  amount  of  an  oxidative  estrogen  compound  that  causes 
depurination  of  guanine  when  it  reacts  with  DNA.  So  while  the  model  included  a  damage 
stage,  it  did  not  relate  a  radiation  dose  to  the  amount  of  damage  inflicted.  Rahmanian  et 
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al.  [80]  also  modeled  the  human  BER  mechanism.  They  utilized  dilFerential  equations 
to  quantify  the  BER  process  in  a  step-by-step  biochemical  kinetic  manner.  In  addition, 
these  researchers  assumed  that  the  total  concentration  of  the  repair  enzymes  was  conserved 
throughout  the  repair  mechanism  in  order  to  normalize  the  differential  equations. 

This  research  will  model  the  damage  and  repair  of  a  spore’s  DNA  following  its 
exposure  to  a  dose  of  ionizing  radiation.  The  primary  DNA  damage  mechanisms  due 
to  ionizing  radiation  are  base  modifications  and  SSBs  which  are  both  repaired  via  the 
NER  process.  Politi  et  al.  [78]  were  among  the  first  to  mathematically  model  the  human 
NER  pathway.  Their  model  considered  damage  via  UV  radiation  and  utilized  a  system 
of  differential  equations  to  sequentially  model  a  six-step  human  NER  pathway.  The 
model  takes  into  account  the  assembly  of  the  necessary  enzymes  and  gives  association 
rates  for  each  enzyme  at  each  step.  The  model  considers  the  assembly  of  each  enzyme 
independently  of  the  other  enzymes  in  the  system  and  regards  the  rate  at  which  each  binds 
to  the  lesion  site  to  be  the  only  important  kinetic  rate  consideration. 

Kesseler  et  al.  [47]  proposed  a  mathematical  model  which  was  also  based  on  the 
association  of  the  proteins  to  the  damage  sites.  It  incorporated  kinetic  “proofreading”  into 
a  human  NER  mathematical  model  and  assumed  that  the  kinetics  in  recruiting  the  correct 
enzymes  to  the  lesion  plays  the  most  important  role.  The  model  integrated  random  order 
assembly  of  repair  factors  and  assumed  any  binding  of  a  repair  enzyme  to  a  damage  site 
will  increase  cooperation  among  repair  factors.  As  in  the  Politi  et  al.  model,  Kesseler  et  al. 
assumed  each  repair  step  takes  the  same  amount  of  time. 

Some  researchers  claim  “repair  mechanisms  can  be  a  double-edge  sword”  because 
simultaneous  repair  can  cause  DSBs  to  occur  when  two  SSBs  lie  relatively  close  together 
on  the  DNA  strand  [46].  Karschau  et  al.  [46]  contend  that  there  exists  a  critical  threshold  of 
the  number  of  repair  enzymes  which  if  larger  than  this  threshold,  the  survival  of  a  bacterial 
cell  is  independent  of  the  number  of  repair  enzymes.  They  relate  cell  death  to  the  rate  of 
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DSB  creation  during  repair  instead  of  direct  or  indirect  DNA  damage.  The  DSB  occurs 
when  one  SSB  is  in  the  process  of  being  repaired  via  NER  and  repair  begins  on  another 
relatively  close  SSB .  Richard  et  al.  [83]  agree  with  this  claim  and  argue  there  is  an  optimum 
level  of  repair  enzymes  in  each  cell  which  is  “explained  as  the  best  trade-off  between  fast 
repair  and  a  low  probability  of  causing”  DSBs  [83].  Mathematical  models  based  on  this 
theory  were  presented  however  neither  model  accounted  for  the  possibility  of  damage  to 
both  the  cell’s  DNA  and  its  repair  enzymes. 

The  amount  of  proteins  available  for  DNA  repair  is  critical  to  the  survival  of  the  spore 
because  a  damaged  spore  will  only  be  able  to  utilize  the  repair  proteins  that  survived  the 
damage  mechanism.  Specifically,  Minton  [60]  showed  that  UvrA  exposed  to  ionizing 
radiation  resulted  in  degraded  DNA  repair.  This  author  introduced  the  concept  of  initial 
protein  ‘fitness’  applied  to  a  population  of  bacteria  and  claimed  that  if  the  level  of  protein 
content  was  degraded  below  a  critical  threshold  level,  the  spore  would  not  have  the  ability 
to  repair  its  DNA  and  thus  would  not  survive  [48].  Hurst  [41]  incorporated  this  idea  of 
initial  protein  ‘fitness’  into  the  Sokhansanj  et  al.  sequential  BER  model  by  assuming  the 
enzymes  experienced  the  same  level  of  damage  as  DNA. 

4.6  Mathematical  Model 

The  repair  model  presented  in  this  research  will  incorporate  the  DNA  damage  model 
described  in  Section  3.8  which  predicts  a  level  of  damage  based  on  a  certain  dose  of 
radiation  thus  coupling  dose,  damage,  and  repair  of  the  Bacillus  spore.  Since  enzymes 
are  necessary  to  repair  DNA,  a  damaged  spore  will  only  be  able  to  germinate  with  the 
enzymes  that  survived  the  radiation  exposure.  DNA  must  be  repaired  before  a  viable  cell  is 
produced  during  outgrowth  and  there  must  be  a  suitable  amount  of  repair  enzymes  available 
to  accomplish  this.  Thus,  the  model  presented  here  will  be  based  on  the  amount  of  enzymes 
available  for  repair  given  that  the  enzymes  were  also  degraded  by  the  ROS  as  described  in 
Section  3.8.3.  In  addition,  this  model  will  include  the  concept  of  an  initial  enzyme  ‘fitness’ 
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density  function  which  was  introduced  by  this  author  and  incorporated  within  the  Hurst 
model. 

Recall  from  Section  4.3  that  there  are  six  enzymes  involved  in  the  NER  process;  UvrA, 
UvrB,  UvrC,  helicase  II,  polymerase  I,  and  ligase.  UvrA  dimerizes  prior  to  the  initiation 
of  the  NER  process.  UvrA2  and  UvrB  form  the  UvrA2Bi  complex  before  binding  to  the 
DNA  damage  site.  Then  UvrA2  conducts  damage  recognition  and  signals  the  damage  to 
UvrB  which  conducts  damage  verification  [24].  Some  researchers  believe  the  complex 
contains  two  UvrB  proteins  which  allows  scanning  of  both  sides  of  the  DNA  to  verify  which 
backbone  contains  the  site  of  the  DNA  damage  [104].  This  NER  method  requires  UvrA2 
and  two  UvrB  enzymes  to  form  a  UvrA2B2  complex  prior  to  binding  to  the  damaged  DNA. 
In  both  cases,  “once  UvrB  has  detected  the  lesion  the  UvrA  subunits  are  released  from  the 
complex”  [37].  Next,  UvrC  attaches  to  the  UvrB-DNA  complex  in  order  to  nick  the  DNA 
strand  on  either  side  of  the  damage  site.  Then,  helicase  II  associates  with  the  damage  site 
and  removes  the  nicked  portion  of  the  DNA  strand  while  also  releasing  the  UvrC  enzyme 
[75].  Polymerase  I  dissociates  both  UvrB  and  helicase  II  from  the  damage  site  and  fills  in 
the  gapped  DNA  with  a  corrected  strand  [15,  69].  Einally,  ligase  releases  the  polymerase  I 
enzyme  and  seals  the  DNA  strand  together. 

This  section  will  offer  a  mathematical  model  for  both  the  UvrA2Bi  and  UvrA2B2  NER 
pathways  since  researchers  do  not  agree  as  to  which  complex  more  accurately  reflects 
the  biological  mechanism  [71].  In  addition,  a  model  which  incorporates  simultaneous 
hydrolysis  damage  with  the  NER  process  will  be  presented.  Einally,  numerical  solutions  to 
these  models  will  demonstrate  the  difference  between  the  two  NER  pathways  and  the  effect 
of  continual  hydrolysis  damage. 

4.6.1  NER  with  UvrA2Bi  Complex. 

The  first  mathematical  model  presented  will  consider  the  NER  pathway  which  utilizes 
the  UvrA2Bi  complex.  The  process  by  which  enzymes  attach  to  and  dissociate  from  the 
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Figure  4.2:  NER  Enzymes  and  Products 


where 


A  -  UvrA 
B  =  UvrB 
C  =  UvrC 
D  =  damaged  DNA 


H  =  helicase  II 
P  =  polymerase  I 
L  -  ligase 
Dr  =  repaired  DNA 


and  the  UvrA2Bi  complex  which  is  fomied  prior  to  its  association  with  the  DNA  lesion  is 
denoted  by  A2B.  Since  UvrA  dimerizes  prior  to  the  start  of  the  NER  process,  this  model 
assumes  that  UvrA  is  released  from  the  DNA  damage  site  as  a  dimer.  This  model  also 
assumes  that  all  helicase  II,  polymerase  I,  and  ligase  enzymes  in  the  bacterial  spore  are 
available  for  the  NER  mechanism.  In  reality,  some  of  these  enzymes  will  participate  in 
other  cellular  functions  such  as  BER. 

This  process  is  described  by  the  following  chemical  reactions: 

A2  +  B  >  A2B  (4. 1 ) 

A2B  +  A2BD  (4.2) 
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A2BD^A2  +  BD  (4.3) 

BD  +  C^CBD  (4.4) 

CBD  +  H^  HBD  +  C  (4.5) 

HBD  +  P^  PD  +  B  +  H  (4.6) 

PD  +  L^  LD  +  P  (4.7) 

LD^Dr  +  L  (4.8) 


where  ki  is  the  reaction  rate  constant  for  each  reaction.  The  chemical  reactions  in  Equations 
4.1  through  4.8  and  enzyme  kinetics  allow  the  NER  mechanism  to  be  modeled  by  the 
following  differential  equations: 


d[A2] 

dt 

d[^] 

dt 

d[C] 

dt 

d[//] 

dt 

d[P] 

dt 

d[L] 

dt 

d[D] 

dt 

d[A2B] 

dt 

d[A2BD] 

dt 

d[fiD] 

dt 

d[CBD] 

dt 

d[HBD] 

dt 

d[ED] 

dt 


=  k3[A2BD]  +  k-dAiB]  -  kdA2][B] 

(4.9) 

=  kdHBD][P]  +  k-dA2B]  -  kdA2][B] 

(4.10) 

=  k5[CBD][H]  -  k4[BD][C] 

(4.11) 

=  kdHBD][P]-k5[CBD][H] 

(4.12) 

=  k7[PD][L]  -  kdHBD][P] 

(4.13) 

=  kdLD]  -  k2[PD][L] 

(4.14) 

=  -k2[A2B][D] 

(4.15) 

=  kdA2][B]  -  k-dA2B]  -  k2[A2B][D] 

(4.16) 

=  k2[A2B][D]  -  k3[A2BD] 

(4.17) 

=  k3[A2BD]  -  k4[BD][C] 

(4.18) 

=  k4[BD][C]  -  k5[CBD][H] 

(4.19) 

=  k5[CBD][H]  -  kdHBD][P] 

(4.20) 

=  ke[HBD][P]-kdPD][L] 

(4.21) 
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(4.22) 


^^=kAPD]{L\-ULD] 
at 

^  =  -t8[LD].  (4.23) 

at 

These  15  differential  equations  eharaeterize  the  proeess  of  repair.  However,  the  system 
size  ean  be  redueed  by  reeognizing  some  eonserved  quantities.  For  example,  adding 
Equations  4.11  and  4.19  gives 

d[C]  ^  djCBD]  ^  ^ 

At  dt 

whieh  implies  [C]  +  [CBD]  is  a  eonstant.  Therefore,  assuming  no  further  damage  is  done 
to  the  enzyme,  the  amount  of  the  UvrC  enzyme  is  always  conserved  whether  the  enzyme  is 
a  part  of  a  complex  or  it  exists  as  a  free  enzyme.  Initially,  at  t  =  0,  the  concentration  of  the 
complex  CBD  is  zero.  Therefore,  let 

[C]  +  [CBD]  =  [C]o  (4.24) 

where  [C]o  represents  the  initial  concentration  of  C.  The  same  enzyme  conservation 
principle  applies  to  helicase  II,  polymerase  I,  and  ligase  since 

d[H]  ^  d[HBD] 

dt  dt 

d[F]  ^  d[FD] 

dt  dt 

d[L]  ^  d[LD] 

dt  dt 

thus  requiring 


[//]  +  [HBD]  =  [H]o 

(4.25) 

[E]  +  [PD]  =  [E]o 

(4.26) 

[L]  +  [LD]  =  [L]o 

(4.27) 

where  [//]o,  [P]o,  and  [L]o  represent  the  initial  concentrations  of  helicase  II,  polymerase  I, 
and  ligase,  respectively.  Substituting  Equation  4.24  into  Equation  4.11  gives 

^  =  ^5  ([C]o  -  [C])  [H]  -  UBDnCl 
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Substituting  Equations  4.24  and  4.25  into  Equation  4.12  yields 

^  =  ^6  ([H]o  -  [//])  [P]  -  ks  ([C]o  -  [C])  [//]. 
at 

Substituting  Equations  4.25  and  4.26  into  Equation  4.13  produces 

^  =  ^7  ([P]o  -  [Z’])  [L]  -  h  i[H]o  -  [//])  [PI 
at 

Substituting  Equations  4.26  and  4.27  into  Equation  4.14  results  in 

^=ks  ([L]o  -  [L])  -  kj  ([E]o  -  [P])  [L]. 
at 

Einally,  substituting  Equation  4.27  into  Equation  4.23  and  Equation  4.25  into  Equation  4. 10 
gives 


^=-t8([L]o-[L]) 

at 

^=ke  ([H]o  -  [H])  [E]  +  k-i[A2B]  -  kdA2][B]. 


These  substitutions  have  eliminated  the  CBD,  HBD,  PD,  and  LD  terms.  Thus  Equations 


4.9  through  4.23  are  reduced  to  the  following: 


d[A2] 

dt 

d[^] 

dt 

d[C] 

dt 

d[H] 

dt 

d[P] 

dt 

d[L] 

dt 

d[D] 

dt 

d[A25] 

dt 

d[A2BD] 

dt 


k3[A2BD]  +  k-dA2B]  -  kdA2][B]  (4.28) 

ke  ([^]o  -  [H])  [E]  +  k.dA2B]  -  kdA2][B]  (4.29) 

k5  ([C]o  -  [C])  [//]  -  k4[BD][C]  (4.30) 

ke  ([//]o  -  [H])  [E]  -  k5  ([C]o  -  [C])  [//]  (4.31) 

ki  ([E]o  -  [P])  [Z.]  -  ke  {[H]o  -  [H])  [E]  (4.32) 

ks  ([L]o  -  [L])  -  kj  ([E]o  -  [P])  [L]  (4.33) 

-k2[A2B][D]  (4.34) 

kdA2][B]  -  k.dA2B]  -  k2[A2B][D]  (4.35) 

k2[A2B][D]  -  k3[A2BD  (4.36) 
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(4.37) 


^^=k^{A2BD]-UBD]{C] 
at 

^^=h([L]o-[L]).  (4.38) 

at 

Note  that  the  amount  of  the  UvrA  enzyme  dimer  is  also  conserved  since  adding  Equations 
4.28,  4.35,  and  4.36  yields 

^  ([A2]  +  [A2B]  +  [A2BD])  =  0. 
at 

Define  [^2)0  to  be  the  initial  concentration  of  A2,  then 

[A2]  +  [A2B]  +  [A2BD]  =  [A2]o 


implies  that 


[A2BD]  =  [A2]o  -  [A2]  -  [A2B]. 


(4.39) 


Similarly,  let  [5]o  be  the  initial  concentration  of  the  UvrB  enzyme.  Adding  Equations  4.10, 
4.16,  4.17,  4.18,  4.19,  and  4.20  gives 

-  ([B]  +  [A2B]  +  [A2BD]  +  [BD]  +  [CBD]  +  [HBD])  =  0 
at 


thus 


[B]  +  [A2B]  +  [A2BD]  +  [BD]  +  [CBD]  +  [HBD]  =  [B]o.  (4.40) 


Substituting  Equations  4.39,  4.24,  and  4.25  into  Equation  4.40  produces 


[BD]  =  ([B]o  -  [B])  -  ([A2]o  -  [A2])  -  ([C]o  -  [C])  -  ([H]o  -  [//]) .  (4.41) 


Equations  4.39  and  4.41  allow  elimination  of  the  A2BD  and  BD  terms.  Therefore, 
Equations  4.28  through  4.38  are  further  reduced  to  the  following  equations, 

^  =h  ([A2]o  -  [A2]  -  [A2B])  +  k-i  [A2B]  -  h  [A2]  [B]  (4.42) 

at 

^  =k(,  ([//]o  -  [H])  [B]  +  k-dA2B]  -  ki[A2][B]  (4.43) 
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^  =k5([C]o-[C])[H] 


-  k,  {([5]o  -  [B])  -  ([A2]o  -  [A2])  -  ([C]o  -  [C])  -  ([H]o  -  [//])}  [C]  (4.44) 

^  =k,  ([H]o  -  [H])  [P]  -  ks  ([C]o  -  [C])  [//]  (4.45) 

at 

^  =k7  ([P]o  -  [P])  [Z^]  -  ke  ([H]o  -  [//])  [P]  (4.46) 

at 

^  =k,  ([L]o  -  [L])  -  kj  ([P]o  -  [P])  [L]  (4.47) 

at 

^=-k2[A2B][D]  (4.48) 

at 

=kdA2][B]  -  k_AA2B]  -  k2[A2B][D]  (4.49) 

at 

^^=k,([L]o-[L]).  (4.50) 

at 

These  equations  are  combined  with  the  initial  conditions,  given  by 

[A2](0)  =[A2]o  (4.51) 

[5](0)  =[B]o  (4.52) 

[C](0)  =[C]o  (4.53) 

[Hm  =[H]o  (4.54) 

[Pm  =[P]o  (4.55) 

[L](0)  =[L]o  (4.56) 

[Dm  =[D]o  (4.57) 

[DM  =0  (4.58) 

[A25](0)  =0,  (4.59) 

to  characterize  the  evolution  of  the  UvrA2Bi  NER  process. 

The  steady  state  solution  of  this  system  (located  in  Appendix  A)  demonstrates  that 
without  additional  damage,  all  of  the  damaged  DNA  will  eventually  be  repaired.  Also 
in  the  long  term,  the  concentrations  of  UvrC,  helicase  II,  polymerase  I,  and  ligase  will 
approach  their  initial  concentrations  as  they  will  no  longer  be  needed  for  the  repair  process. 
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In  addition,  the  concentrations  of  UvrA,  UvrB,  and  the  UvrA2Bi  complex  will  exist  in  their 
equilibrium  state  given  by  the  equation  k-i[A2E\  =  ki[A2][B]  [85]. 

4.6.2  NER  with  UvrA2B2  Complex. 

In  some  cases  of  NER,  the  UvrA  dimer  forms  a  complex  with  two  UvrB  molecules  as 
shown  by  Figure  4.3  [37].  The  UvrA2B2  complex  allows  UvrB  to  scan  both  DNA  strands 
simultaneously  in  order  to  verify  the  damage  detected  by  UvrA.  UvrB  does  not  form  a 
dimer  during  this  process.  Instead,  the  two  UvrB  molecules  attach  to  either  side  of  the 
UvrA  dimer  [72].  Since  UvrB  does  not  form  a  dimer,  the  UvrB  molecules  are  released  as 
two  monomers  by  the  polymerase  I  enzyme. 


Figure  4.3:  NER  with  UvrA2B2  Complex 


This  variation  of  the  NER  pathway  is  defined  by  the  following  chemical  reactions: 


A2  "1"  B 

h 

^  A2B 

k-i 

(4.60) 

A2B  +  B 

h 

,  ■  ^  A2B2 

k-i 

(4.61) 

A2B2  H"  D 

ki 

-A  A2B2D 

(4.62) 

A2B2D 

— >  A2  +  B2D 

(4.63) 
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B2D  +  C- 

CB2D 

(4.64) 

CB2D  +  H- 

HB2D  +  C 

(4.65) 

HB2D  +  P- 

PD +  2B  +  H 

(4.66) 

PD  +  L- 

h 

LD  +  P 

(4.67) 

LD- 

Df  +  L. 

(4.68) 

These  chemical  reactions  are  also  modeled  by  a  system  of  differential  equations. 


^  =  kAA2B2D]  +  k.AA2B]  -  kdA2][B] 
at 

^  =  2h[HB2D][P]  +  k_dA2B]  +  k.dA2B2]  -  h[A2][B]  -  h[A2B][B] 
at 

^  =  k5[CB2D][H]  -  k4[B2D][C] 
at 

^  =  h[HB2D][P]  -  k5[CB2D][H] 
at 

^  =  k2[PD][L]  -  h[HB2D][P] 
at 

d[L] 

-^=ks[LD]-k2[PD][L] 

at 

d[D] 

=  -k2[A2B2][D] 
dt 

d[A25] 

-4^  =  kdA2][B]  +  k-dA2B2]  -  k.dA2B]  -  ki[A2B][B] 
dt 

d[A2B2]  ~ 

=  kdA2B][B]  -  k.dA2B2]  -  k2[A2B2][D] 


d[A2B2D] 

dt 


=  k2[A2B2][D]  -  k3[A2B2D] 


d[B2D] 

dt 


=  kAA2B2D]  -  k4[B2D][C] 


d[CB2D] 

dt 

d[HB2D] 

dt 

d[PD] 

dt 

d[LD] 

dt 

d[T>.] 

dt 


=  k4[B2D][C]-k5[CB2D][H] 

=  k5[CB2D][H]  -  kdHB2D][P] 
=  kdHB2D][P]  -  kdPD][L] 

=  kdPD][L]  -  kdLD] 

=  kdLD]. 


(4.69) 

(4.70) 

(4.71) 

(4.72) 

(4.73) 

(4.74) 

(4.75) 

(4.76) 

(4.77) 

(4.78) 

(4.79) 

(4.80) 

(4.81) 

(4.82) 

(4.83) 

(4.84) 
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The  amount  of  each  enzyme  in  the  NER  process  that  utilizes  the  UvrA2B2  complex  is  still 
conserved  whether  the  enzyme  is  part  of  a  compound  or  it  exists  as  a  free  enzyme.  Note 
that  the  stochiometry  of  the  UvrB  enzyme  in  this  system  leads  to  a  different  conservation 
principle  for  the  UvrB  molecule: 

d 

—  ([B]  +  [A2B]  +  2[A2B2]  +  2[A2B2D]  +  2[B2D]  +  2[CB2D]  +  2[HB2D])  =  0. 

dt 

Therefore  let 


[B]  +  [A2B]  +  2[A2B2]  +  2[A2B2D]  +  2[B2D]  +  2[CB2^]  +  2[HB2D]  =  [B]o 


Using  this  conservation  equation  in  addition  to  the  conservation  equations  previously 
developed  in  Equations  4.24  through  4.27  and  4.39,  Equations  4.69  through  4.84  are 
reduced  to  the  following  ten  equations, 


^  =  h  ([A2]o  -  [A2]  -  [A2B]  -  [A2B2])  +  k-dA2B]  -  h[A2][B] 
dt 

^  =  2h  ([//]o  -  [H])  [B]  +  k.dA2B]  +  Li[A2B2]  -  kdA2][B]  -  h[A2B][B] 
dt 

^  =  ^5  ([c]o  -  [c])  [//]  -  k4 ([B]o  -  m 

-  ([A2]o  -  [A2])  +  ^[A2B]  -  ([C]o  -  [C])  -  ([//]o  -  [//])|  [C] 

^  =  ^6  ([//]o  -  [H])  [B]  -  k5  ([C]o  -  [C])  [H] 
dt 

^=k2  ([B]o  -  [B])  [L]  -  h  ([//]o  -  [H])  [B] 
dt 

^  =  ^8  ([Z.]o  -  [L])  -  kj  ([B]o  -  [B])  [L] 
dt 

d[D] 

=  -k2[A2B2][D] 
dt 

d[A2B] 

-4^  =  kdA2][B]  +  k.dA2B2]  -  k-i[A2B]  -  kdA2B][B] 

dt 

d[A2B2]  ~ 

-44^  =  kdA2B][B]  -  k.dA2B2]  -  k2[A2B2][D] 

dt 

d[Dd 

-^=k,i[L]o-[L]). 


(4.85) 

(4.86) 

(4.87) 

(4.88) 

(4.89) 

(4.90) 

(4.91) 

(4.92) 

(4.93) 

(4.94) 
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Again,  we  combine  these  equations  with  the  initial  conditions  given  in  Equations  4.51 
through  4.58  with  the  addition  of  the  initial  condition  [A252](0)  =  0,  to  characterize  the 
evolution  of  the  UvrA2B2  repair  process. 

4.6.3  Characteristic  Scaling. 

To  convert  the  variables  into  dimensionless  variables,  certain  characteristic  scalings 
are  introduced.  This  allows  all  calculations  to  be  independent  of  particular  units  chosen.  In 
order  to  non-dimensionalize  enzyme  concentration  and  time,  let 


[A2](0  =  [A2]o^i(t) 

(4.95) 

[B](0  =  [B]o6(r) 

(4.96) 

[C](0  =  [C]o^3(t) 

(4.97) 

[Hm  =  [H]oUr) 

(4.98) 

=  [PUsir) 

(4.99) 

[im  =  [LUeir) 

(4.100) 

[Dm  =  [DT]Ur) 

(4.101) 

[D^m  =  [Drmir) 

(4.102) 

where  r  =  k^t  so  t  is  dimensionless  time  and  Dj  is  the  total  concentration  of  DNA  in  the 
healthy,  undamaged  spore.  These  definitions  were  chosen  because  the  initial  concentrations 
of  each  enzyme  are  the  maximum  value  each  free  enzyme  concentration  can  be.  Also,  the 
maximum  concentration  that  both  damaged  and  repaired  DNA  can  have  is  the  concentration 
of  the  total  DNA  in  the  spore.  The  time  scale,  k^,  was  chosen  because  this  rate  parameter 
has  units  of  sec“^  Finally,  let 

[A2Bm  =  [A2]o6(t)  (4.103) 

as  the  value  of  [A2]o  represents  the  maximum  concentration  that  A2B  can  be.  Note  that 

=  t:  ([^2]o^i(t))  =  \A2\o-r~^  = 
at  at  at  dr  dr 
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In  order  to  completely  non-dimensionalize  the  system,  define  dimensionless  parameters 


= 


h 


h  = 


h[B]c 


h  = 


kiiPh 

k-i 


a\ 


in 

[B]c 


ai 


me 


m]o 


0^3 


[A2]o 


a4 


A4  = 

[H]o 


ksiCh 
h  ’ 


^5  = 


ki[L\ 


[C]o’ 
and  jS3 


mo’ 

^  o  -  ^-1 

4.6.3. 1  Non-dimesionalized  NER  UvrA2Bi  Pathway. 


[Pk 

me 


CK6 


ki‘ 


[A2]o 

[Dr]’ 

(4.104) 


Utilizing  Equations  4.95  through  4.103  and  the  dimensionless  parameters,  Equations 
4.42  through  4.50  become: 


dr 

d^2 

dr 

d6 

dr 


=  l-^l+^7(;S2-l)-d2^l6 
=Aiai  (1  -  ^4)^5  +  cr3^7^2  -  A2as^i^2 
=d4«4(l  -  ^3)^4 


-  42^3^3 


1  -  ^2  -  0^3  (1  -  ^1) - ^  (1  -  ^3) - ^  (1  -  ^4) 

a2Q;4  Q'2 


d^4 

dr 

dr 

d^6 

dr 

d^d 

dr 

d^i 

dr 


=4ia2(l-^4)^5 -44(1 -6)^4 


=45(1  -  ^5)^6 -4i(l -^4)^5 

=A(l-^6)-45a5(l-6)^6 

=  -  43Q;6^7^rf^2 

=42^1^2  -  ^1^2  (1  +  43^af) 


dr 


A«6Q'i 

0:30^5 


(1-^6). 


(4.105) 

(4.106) 

(4.107) 

(4.108) 

(4.109) 

(4.110) 

(4.111) 

(4.112) 

(4.113) 


This  numerical  solution  to  this  system  of  non-dimensional  differential  equations  will  be 
presented  in  Section  4.6.5. 
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4.63.2  Non-dimesionalized  NER  UvrA2B2  Pathway. 


In  order  to  create  a  dimensionless  system  for  the  alternate  NER  pathway  that  utilizes 
the  UvrA2B2  complex,  additional  dimensionless  parameters  must  be  defined.  Let 


kdH]o 

h 


ki[B]o 


and 


^2  = 


k-i 

h  ' 


In  addition,  denote 


{A2B2m  =  [A2]o^8(t) 


(4.114) 


in  order  to  non-dimensionalize  concentration  and  time.  With  Equations  4.95  through  4.102, 
Equation  4.114,  and  the  dimensionless  parameters.  Equations  4.85  through  4.94  become: 


-l\ai  (1  -  ^4)  ^5  +  Q'3  (^7^2  +  ^8^2)  -  ^3^2  (^2^1  + 


=/l4«4(l  -  ^3)^4 


^=1-^1  (1+42^2)  -^8  -^7(1  ->02) 
dr 

d^2 

dr 
dr 

d^4 
dr 

=45(1  -  ^5)^6  -  di  (1  -  ^4)  ^5 
dr 


1 


(l-^2)-a3ll-^i-k7 


^2/33^3 


=4ia2  (1  -  ^4)  ^5  -  44  (1  -  ^3)  ^4 


^(1-^3)--(1-^4) 
a2Q^4  ®2 


d^6 

dr 


=A(l-^6)-d5a5(l-^5)^6 


-  -  43Cr6^8^rf/32 

f  =42^16  +^8.02  -6  (>02  + te) 

^=fe6-^8(>02  +  43^,;02) 

/3iaeai 

-r-  = - -«)• 

dr  0:30:5 


(4.115) 

(4.116) 

(4.117) 

(4.118) 

(4.119) 

(4.120) 

(4.121) 

(4.122) 

(4.123) 

(4.124) 


This  system  is  also  solved  numerically  and  portrayed  in  Section  4.6.5. 


4.6.4  Hydrolysis  Damage. 

This  repair  model  assumes  that  all  DNA  damage  is  caused  by  radiation  and  occurs 
prior  to  germination.  However,  if  a  spore  receives  the  nutrients  it  needs  to  germinate. 
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then  the  spore  core  will  take  in  water  which  can  produce  hydrolysis  damage.  Hydrolysis 
is  a  reaction  of  a  molecule  with  water.  Water  can  react  with  polymers;  such  as  DNA 
and  enzymes,  resulting  in  depolymerization  (breaking  a  large  strand  of  polymers  into  two 
smaller  ones)  or  removal  of  a  side  group  of  the  polymer  [58]. 

In  addition,  water  can  react  with  DNA  causing  base  modifications  such  as  deamination 
and  depurination  as  discussed  in  Section  3.6.1.  These  base  modifications  create  a  structural 
change  in  the  DNA  backbone  which  prevents  the  DNA  from  replicating  properly  thus 
causing  spore  death.  This  section  will  introduce  a  repair  model  which  also  includes 
hydrolysis  damage  to  the  DNA  only. 

Define  the  undamaged  DNA  in  the  cell  or  “good  DNA”  as  Dg.  Because  damaged 
DNA  is  compounded  with  the  repair  enzymes  during  the  NER  process,  the  total  DNA  in 
the  spore  at  any  time,  t,  is  not  equal  to  the  sum  of  the  damaged  DNA  and  the  viable  DNA. 
Rather,  the  total  DNA  is  a  combination  of  the  undamaged  DNA,  damaged  DNA  that  has 
not  begun  to  be  repaired,  and  DNA  that  is  in  the  process  of  being  repaired,  i.e. 

[Dt]  =  [DgM  +  [D](0  +  [DA(t) 

where  [Dc](0  represents  damaged  DNA  that  is  compounded  with  repair  enzymes.  That  is, 

[D,](0  =  [AiBOm  +  [HBD](t)  +  [CBD](t)  +  [BD](t)  +  [PD](t)  +  [LD]it). 

Note  that,  [Dg](0)  =  [Dj]  -  [D]o  because  [Dc](0)  =  0.  With  continual  hydrolysis  damage 
to  the  DNA,  the  NER  model  in  Equations  4.42  through  4.50  require  the  addition  of  a 

differential  equation  for  [D^], 

d[DJ 

=  h  ([^]o  -  [L])  -  k*[Dg],  (4.125) 

dt 

where  k*  is  the  reaction  rate  constant  that  represents  the  averaged  rate  for  all  possible 
hydrolysis  reactions  in  the  spore.  The  only  other  change  is  the  addition  of  one  term  to 
Equation  4.48  which  becomes 

^  =  -k2[A2B][D]  +  k*[Dg].  (4.126) 
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These  same  two  alterations  apply  to  both  the  NER  model  with  the  UvrA2Bi  complex  and 
the  model  with  the  UvrA2B2  complex  in  Equations  4.85  through  4.94.  The  value  k*  follows 
the  Arrhenius  equation  which  models  the  rate  of  chemical  reactions  [58]: 

k*{T)  =  Ae-^{H20] 

where  [H2O]  is  the  concentration  of  water.  Temperature,  T ,  is  assumed  to  be  room 
temperature  (i.e.  300  °K)  and  the  gas  constant,  R,  is  1.986x10“^  kcal/mol°K.  The  encounter 
frequency,  A  and  the  activation  energy,  E,  are  dependent  upon  the  moisture  content  of  the 
spore.  Values  of  k*  for  various  levels  of  moisture  content  are  displayed  in  Table  4.1. 


Table  4.1:  Parameters  for  k*  [108] 


* 

1 

Moisture  Content  (%) 

E  (kcal/mol) 

A*[H20]  (s-^) 

3.77x10-'^ 

5 

26.1 

4.0x10^0 

1.38x10-* 

9 

24.5 

1.0x10^° 

1.95x10-5 

19 

12.3 

1.8x10^ 

8.2x10-5 

30 

10 

1.6x105 

Similarly,  characteristic  scaling  is  introduced  in  order  to  create  a  dimensionless  system 
that  represents  the  continual  hydrolysis  damage.  With 

[D,](t)  =  [Dr]4(T)  and  J84  = 

^3 


along  with  the  dimensionless  parameters  and  variables  utilized  in  the  original  repair  model 
(Equations  4.95  through  4.104),  Equations  4.125  and  4.126  become 


d^g  pia^ai 


(1-^6) -Mr 


dr  a^as 

The  mathematical  NER  model  that  includes  hydrolysis  damage  is  also  implemented 
and  solved  numerically.  The  results  will  be  demonstrated  in  the  next  section. 
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4.6.5  Results. 


The  dimensionless  repair  models  are  numerically  evaluated  via  a  moderately  stiff 
MATLAB  solver  using  rate  kinetics  and  enzyme  concentrations  found  in  the  literature. 
A  stiff  solver  is  required  for  this  system  of  differential  equations  because  the  time  scales 
differ  widely  and  “the  solution  being  sought  varies  slowly,  but  there  are  nearby  solutions 
that  vary  rapidly”  [63].  The  rate  parameters  are  reported  in  Table  4.2.  Recall  from  Section 
3.8.2  that  the  E.  coli  cell’s  DNA  concentration  prior  to  radiolysis  is  2.87  x  10^  nM  therefore 
this  value  is  applied  to  {Dj].  The  same  section  proved  that  10  percent  of  the  initial  DNA  is 
damaged  therefore  [D]o  =  .1  *  [Dj-].  Enzyme  concentrations  prior  to  radiation  damage  are 
found  in  Table  3.6.  As  the  enzymes  were  also  reduced  to  90  percent  of  their  initial  activity 
level  via  radiation  damage,  the  amount  of  each  enzyme  available  for  the  repair  process  after 
a  10  percent  reduction  is  listed  in  Table  4.3. 


Table  4.2:  Rate  Parameters 


Parameter 

Value 

Notes 

ki 

1.1x10-’  nM-^s-i 

Geobacillus  stearothermophilus  NER  [72] 

k_i 

8.5x10-5  s-i 

NERofE.  coli  [70] 

k2 

6x10-5  nM-^s-^ 

UV-irradiated  E.  coli  [69] 

ks 

3x10-5  s-^ 

UV-irradiated  E.  coli  [100] 

k4 

8.3x10-*  nM-^s-i 

E.  coli  NER  [109] 

k5 

0.12nM-is-^ 

UvrD  in  E.  coli',  assumes  UvrD  dimer  [55] 

ke 

6.2x10-^^  nM-^s"^ 

E.  coli  BER  of  DNA-methylation  damage  [101] 

ky 

8.5x10-“  nM-^s-^ 

E.  coli  BER  of  DNA-methylation  damage  [101] 

ks 

0.075  s-^ 

UV-irradiated  E.  coli  [52] 
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Table  4.3:  Data  for  E.  coli  after  Damage 


Molecule 

Concentration 

[A2]o 

17.2  nM 

[B]o 

336.2  nM 

[C]o 

14.9  nM 

[//]o 

4109.9  nM 

[P]o 

657.5  nM 

[L]o 

298.9  nM 

[D]o 

2.87  X  10^  nM 

Figure  4.4  displays  the  numerical  results  and  differences  of  the  NER  process  when 
utilizing  the  two  UvrA2B  complexes.  The  horizontal  axis  represents  real  time  in  hours 
and  the  vertical  axis  corresponds  to  the  fraction  of  DNA  repaired  so  that  0.1  represents  all 
damaged  repaired.  The  blue  line  represents  repair  over  time  of  a  DNA  strand  via  NER 
with  the  UvrA2Bi  complex.  This  result  is  consistent  with  conclusions  from  researchers  as 
Moeller  et  al.  [61]  found  that  81  percent  of  SP  and  CPDs  in  a  B.s.  spore  were  repaired 
within  one  hour  (recall  that  SP  and  CPDs  are  repaired  via  NER).  The  dotted  red  line 
represents  the  NER  process  that  utilizes  the  UvrA2B2  complex  and  assumes  that  ki  and 
k\  are  equal  in  value.  Assuming  the  same  amount  of  initial  DNA  damage,  repair  with  the 
UvrA2B2  complex  is  slower  than  that  with  the  UvrA2Bi  complex  because  of  the  increased 
requirement  for  the  UvrB  enzyme.  However,  ki  is  most  likely  smaller  in  value  than  ki  so 
results  of  the  repair  process  with  the  UvrA2B2  complex  and  a  ki  value  that  is  one  order  of 
magnitude  smaller  than  that  of  k\  are  plotted  via  the  black  dashed  line.  The  decrease  in  the 
k\  rate  causes  NER  with  the  UvrA2B2  complex  to  need  an  even  longer  period  of  time  to 
complete  its  DNA  repair  process. 
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UvrA^B  vs  UvrA-B^ 

2  2  2 


Figure  4.5  depicts  the  behavior  of  UvrA,  UvrB,  and  the  UvrA2Bi  complex  during  the 
NER  mechanism.  The  horizontal  axis  corresponds  to  time  in  hours  and  it  represents  the 
same  time  scale  as  the  horizontal  axis  in  Figure  4.4.  The  vertical  axis  gives  the  fraction  of 
enzyme  concentration  available  during  the  NER  process.  The  inset  in  Eigure  4.5  represents 
the  initial  behavior  of  UvrA  and  UvrB  in  blue  and  red  respectively.  The  behavior  reflected 
in  the  inset  graph  occurs  within  the  first  second  of  the  initiation  of  NER.  Note,  that  both 
enzyme  concentrations  decrease  initially  as  the  UvrA2Bi  complex  is  formed.  In  addition, 
UvrA  decreases  at  a  faster  rate  than  UvrB  as  there  is  a  lower  concentration  of  UvrA  in 
the  spore  than  UvrB  (see  Table  4.3).  However,  UvrA  is  released  from  the  damage  site  as 
soon  as  it  passes  the  damage  to  UvrB  for  verification  and  then  becomes  available  to  form 
the  UvrA2Bi  complex  again  very  quickly.  UvrB  is  bound  to  the  damaged  DNA  until  it  is 
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released  by  polymerase  (refer  baek  to  Figure  4.2)  and  it  is  required  to  form  the  UvrA2Bi 
eomplex  as  soon  as  it  is  released.  Therefore,  UvrB  is  utilized  almost  at  its  maximum 
eapacity  during  the  NER  proeess.  The  initial  eoneentration  of  the  UvrA2Bi  complex 
is  assumed  to  be  zero.  The  complex  attaches  to  the  damage  site  almost  immediately 
after  it  is  formed.  Note  that  when  the  majority  of  the  damaged  DNA  is  repaired,  UvrA, 
UvrB,  and  the  UvrA2Bi  complex  approach  their  equilibrium  states  given  by  the  equation 
k-AAiB]  =  ki[A2\m  [85]. 


0  0.5  1  1.5  2  2.5  3 


Time  (hrs) 

Figure  4.5:  UvrA,  UvrB,  and  UvrA2B  Behavior  During  NER 


Helicase  and  polymerase  behavior  during  the  NER  process  is  demonstrated  in  Eigure 
4.6.  The  range  of  time  of  the  horizontal  axis  is  the  same  as  the  previous  two  graphs. 
Notice  that  the  vertical  axis  indicates  almost  45  percent  of  polymerase  is  bound  to  the  DNA 
damage  site  for  repair  while  approximately  8  percent  of  helicase  is  utilized  for  repair.  The 
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difference  is  due  to  the  varying  initial  concentrations  of  each  enzyme  (see  Table  4.3).  The 
inset  of  Figure  4.6  represents  the  initial  behavior  of  helicase  and  polymerase  which  occurs 
within  two  minutes  after  the  start  of  NER.  Recall  that  helicase  is  required  for  NER  prior  to 
polymerase.  Interestingly,  the  maximum  fraction  of  UvrB,  helicase,  and  polymerase  that 
is  bound  to  the  DNA  damage  site  translates  to  an  equivalent  concentration  for  each  of  the 
three  enzymes.  This  indicates  that  the  concentration  of  UvrB  dictates  the  concentrations  of 
helicase  and  polymerase  that  are  bound  to  the  damage  site. 


Helicase  and  Polymerase 


The  amount  of  ligase  bound  to  the  DNA  damage  site  during  NER  is  shown  in  Eigure 
4.7.  The  vertical  axis  corresponds  to  fraction  of  ligase  used  during  the  repair  process.  The 
behavior  of  ligase  over  time  mimics  the  behavior  of  the  helicase  and  polymerase  enzymes. 
However,  only  3.3x10“’  percent  of  the  initial  ligase  concentration  is  bound  to  the  DNA 
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damage  site  beeause  ligase  is  only  attached  for  one  step  of  the  NER  mechanism  and  the  rate 
parameter  for  this  step  is  relatively  fast  compared  to  other  rate  constants  (see  Table  4.2). 
Rahmanian  et  al.  [80]  studied  of  the  accumulation  of  enzymes  to  the  DNA  damage  site 
during  BER  in  human  cells.  Eukaryotic  BER  employs  enzymes  from  both  the  polymerase 
and  ligase  families.  The  qualitative  behavior  of  the  polymerase  and  ligase  enzymes  as 
they  accumulate  to  the  damage  site  in  this  research  is  extremely  similar  to  the  laboratory 
results  found  by  Rahmanian  et  al.  Similar  to  ligase,  only  2.5x10“^  percent  of  UvrC’s  initial 
concentration  is  bound  to  the  damage  site  during  NER.  UvrC  is  also  only  attached  to  the 
DNA  for  one  step  of  the  repair  process. 


Eigure  4.8  depicts  the  result  of  simultaneous  hydrolysis  damage  on  the  repair  model. 
The  graph  portrays  the  effects  of  differing  values  for  k*  on  the  production  of  healthy 
DNA  over  time  given  that  ten  percent  of  the  DNA  is  damaged  prior  to  initiation  of 
repair.  Specifically,  the  vertical  axis  represents  the  percentage  of  undamaged  DNA  and 
the  horizontal  axis  corresponds  to  time  in  hours.  The  solid  black  line  demonstrates  the 
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scenario  with  no  hydrolysis  damage  (i.e.  k*=0).  The  red  line  represents  the  pereentage 
of  DNA  repaired  over  time  given  a  moisture  eontent  of  5  pereent  (refer  to  Table  4. 1  for 
the  eorresponding  k*  values).  Sinee  hydrolysis  damage  is  eontinuously  oeeurring,  DNA  is 
never  100  pereent  repaired.  A  5  pereent  moisture  eontent  will  permit  repair  of  almost  all  of 
the  damaged  DNA  (99.75  pereent  of  the  total  DNA  is  viable).  As  shown  by  the  dashed  blue 
line,  a  moisture  eontent  of  9  pereent  requires  an  even  longer  period  of  time  for  99  pereent 
of  the  DNA  to  be  healthy.  These  results  from  the  hydrolysis  model  are  also  eonsistent 
with  experiments.  Uphoff  et  al.  [101]  found  that  oxidative  DNA  damage  eaused  by  the 
normal  metabolism  of  E.  coli,  required  approximately  two  hours  to  repair.  If  the  rate  of 
hydrolysis  reaetions  overwhelms  the  repair  eapabilities  of  the  spore,  the  eoneentration  of 


Hydrolysis  Damage 


undamaged  DNA  will  deerease  over  time.  This  seenario  is  illustrated  by  the  dotted  blaek 
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line  in  Figure  4.8  whieh  represents  repair  with  a  k*  that  eorresponds  to  approximately  a  10 
pereent  moisture  eontent.  Therefore,  as  k*  inereases,  the  rate  at  whieh  DNA  is  repaired  will 
deerease.  In  addition,  the  pereentage  of  viable  DNA  that  is  repaired  will  deerease.  This 
trend  eontinues  until  the  rate  of  hydrolysis  damage  overwhelms  the  spore’s  repair  eapaeity 
and  the  amount  of  healthy  DNA  eontinuously  deereases  as  shown  by  the  dotted  blaek  line. 

4.7  Summary 

This  ehapter  developed  and  examined  repair  proeesses  for  nueleotide  damage  to  DNA. 
The  mathematieal  model  eaptures  the  relationship  between  six  essential  enzymes  and  the 
repair  of  damaged  DNA.  The  ehemieal  kineties  elearly  show  that  the  loss  of  any  one 
enzyme  stops  the  repair  proeess.  It  follows  that  the  reduetion  of  them  will  slow  repair. 

Two  distinet  pathways  to  repair  requiring  either  UvrA2Bi  or  UvrA2B2  complexes 
were  examined.  Results  indicate  that  the  UvrA2Bi  pathway  was  faster  than  the  UvrA2B2 
pathway.  This  was  likely  due  to  the  added  step  in  formulating  UvrA2B2  as  well  as  the 
reduction  in  availability  of  the  UvrB  molecule. 

The  UvrA2Bi  pathway  was  further  extended  to  include  the  possibility  of  hydrolysis 
damage  during  the  repair  process.  Because  hydrolysis  damage  occurs  simultaneously  with 
repair,  the  viable  DNA  may  never  reach  the  necessary  threshold  for  survival.  The  influence 
of  moisture  content  is  illustrated  in  Figure  4.8.  The  initial  ‘fitness’  of  each  spore  is  reflected 
in  the  state  of  the  essential  enzyme  population  as  well  as  the  DNA.  With  the  addition  of 
hydrolysis  damage,  if  this  initial  enzyme  viability  is  not  sufficient  for  DNA  repair,  the 
population  will  not  survive.  The  probability  of  survival  of  a  spore  population  is  modeled  in 
the  next  chapter. 


80 


V.  Probability  of  Survival 


The  previous  two  ehapters  have  foeused  on  the  damage  and  repair  of  a  single  spore. 
This  ehapter  will  eonsider  the  probability  of  survival  for  a  population  of  spores  assuming 
the  population  has  a  eertain  initial  enzyme  ‘fitness’  as  described  in  Section  3.8.3.  First, 
previous  literature  regarding  biological  survival  curves  will  be  reviewed.  Then,  the 
methodology  for  the  probability  of  survival  model  developed  for  this  research  will  be 
introduced,  and  finally,  results  will  be  presented. 

5.1  Literature  Review 

Past  mathematical  modeling  of  biological  survival  curves  has  centered  on  data-fitting 
of  probability  of  survival  curves  because  there  exists  experimental  data  that  can  be  used  for 
comparison.  For  example,  Brahmi  et  al.  [12]  modeled  the  UV  inactivation  of  bacteria  in 
wastewater  as 

m  _  -Ku 
No 

where 


N(t)  =  number  of  bacteria  at  time,  t 
No  =  number  of  initial  bacteria 
K  =  rate  constant 
I  =  UV  intensity 
t  =  time 

and  instituted  small  deviations  based  on  experimental  results.  This  model  is  also  known  as 
the  single  stage  exponential  decay  equation. 

Prior  to  Brahmi  et  al.’s  work,  Kowalski  et  al.  [49]  used  the  same  model  to  demonstrate 
the  disinfection  of  air  using  UV  irradiation.  These  researchers  adapted  the  model  to 
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consider  two  microbial  populations;  one  that  decays  with  initial  exposure  and  one  that 
is  more  resistant  and  will  survive  for  longer  periods  of  time.  The  model  is 

where  S  (t)  is  the  surviving  fraction  of  the  initial  microbial  population  at  time,  t,  f  is  the 
resistant  fraction  of  the  population,  and  ki  and  ^2  are  rate  constants.  This  equation  is  known 
as  the  two-stage  model.  Kowalski  et  al.  also  modified  their  models  to  account  for  what 
is  known  as  the  quasi-threshold.  The  quasi-threshold  gives  the  shoulder  of  a  biological 
survival  curve  that  represents  the  period  of  time  prior  to  which  a  bacterial  population 
experiences  exponential  decay.  In  other  words,  there  exists  some  accumulation  of  damage 
before  the  onset  of  exponential  decay  which  is  not  demonstrated  until  the  quasi-threshold 
time  period  is  complete.  This  shouldered  survival  curve,  as  shown  in  Figure  5.1  is  also 
known  as  the  sigmoid  survival  curve. 


Dose 


Figure  5.1:  Sigmoid  Survival  Curve 


The  single  stage  and  two  stage  models  are  extremely  similar  to  the  models  developed 
within  the  target  theory  of  cell  survival.  Target  theory  assumes  that  a  cell  contains  a  small 
volume  of  targets  which  can  absorb  energy  to  produce  damage,  otherwise  known  as  a  hit. 
In  target  theory,  the  shoulder  portion  of  the  survival  curve  represents  the  number  of  targets 
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that  had  to  be  inactivated  before  outgrowth  of  a  bacterial  colony  can  be  deterred.  The 
single-hit  model  assumes  that  the  cell  can  only  survive  if  it  sustains  no  hits  and  it  will  die 
if  it  experiences  one  or  more  hits.  The  single-hit  model  is  represented  by 

S(d)  = 

where  S  is  the  survival  probability  per  cell,  d  is  the  dose,  or  hits  per  unit  volume,  and  do  is 
the  dose  required  to  reduce  the  surviving  fraction  to  ^  of  its  initial  value.  This  model  does 
not  account  for  the  quasi-threshold,  thus  it  was  adapted  to  form  a  multitarget-single-hit 
model: 

S{q,n,d)  =  1  -(l  -e-''^)” 

where  each  target  is  assumed  to  have  the  same  probability  of  being  hit  {q)  and  there 
exists  n  targets  in  each  cell.  Note  that  for  n  greater  than  one,  the  multitarget- single¬ 
hit  model  requires  a  zero  slope  when  the  dose  is  equal  to  zero.  Experimental  data  has 
not  replicated  this  phenomenon,  therefore  a  modified  multitarget-single-hit  model  was 
developed  to  include  the  assumption  that  some  fraction  of  cells  irradiated  are  killed  via 
direct  action.  This  leads  to 

5  (^,  n,  J)  =  [  1  -  ( 1  -  e"^"^)”] 

where  q^  is  the  inactivation  coefficient  for  the  cell  killing  that  arises  from  single  hits  and 
qn  is  the  inactivation  coefficient  for  the  previous  muititarget-singie-hit  modei  [4].  Katz  [9] 
ciaimed  the  modified  muititarget-singie-hit  modei  demonstrated  the  ceii’s  inactivation  from 
a  dose  with  two  different  modes  such  as  iow  LET  and  high  LET. 

Target  theory  is  a  reasonabie  approach  to  modeiing  the  survivai  of  prokaryotic  ceiis. 
If  a  dormant  spore  experiences  radiation  damage,  it  wifi  not  begin  its  repair  processes 
untii  it  germinates  and  has  a  vegetative  metaboiism.  This  negates  the  need  to  modei  the 
simuitaneous  repair  of  damage.  However,  modeiing  of  eukaryotic  ceii  survivai  requires 
consideration  of  these  concurrent  procedures.  Therefore  the  iinear-quadratic  moiecuiar 
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model  was  formulated  by  Chadwiek  and  Leenhouts  to  aeeount  for  simultaneous  enzymatie 
repair  [4].  Another  benefit  of  the  linear-quadratie  model  is  that  it  aeeounts  for  the  likelihood 
of  two  SSBs  beeoming  a  DSB  when  they  oeeur  at  a  similar  time  and  in  elose  proximity. 
This  parameter  is  referred  to  as  the  biologieal  effeetiveness  faetor  and  is  denoted  by  p  in 
the  linear-quadratie  model  given  by 

S{d)  = 

where  a  is  a  funetion  of  the  number  of  sites  that  ean  sustain  a  DSB,  the  hit  probability 
eonstant,  the  unrestored  fraetion  of  DSBs,  and  the  fraetion  of  the  dose,  d,  that  ruptures  both 
strands  in  a  single  hit.  The  parameter  jS  is  a  funetion  of  the  unrestored  fraetion  of  DSBs,  the 
effeetiveness  faetor  of  two  SSBs  beeoming  a  DSB,  the  number  of  eritieal  moleeular  bonds 
on  eaeh  DNA  strand,  the  unrestored  fraetion  of  bonds  in  eaeh  strand,  the  hit  probability 
eonstant,  and  the  fraetion  of  the  dose  that  ruptures  one  of  the  strands  in  a  single  hit.  In 
order  to  gain  the  benefit  of  the  effeetiveness  faetor  and  apply  the  linear-quadratie  model 
to  the  survival  of  prokaryotie  eells  that  do  not  allow  simultaneous  repair  of  damage,  the 
unrestored  fraetion  of  breaks  ean  be  set  to  one  [4] . 

Brynjolfsson  [13]  transformed  the  target  theory  models  by  assuming  that  a  mierobe 
in  a  eell  that  is  essential  for  survival  has  a  eomplement,  sueh  as  double-stranded  DNA. 
He  further  assumed  this  essential  mierobe  is  eomposed  of  subunits  and  if  one  subunit  is 
destroyed,  then  its  eomplement  ean  be  used  to  reereate  it.  While  this  model  was  validated 
by  experimental  data,  it  assumes  that  if  both  eomplementary  subunits  are  hit,  the  eell  eannot 
survive.  Seetion  4.4  diseussed  the  ability  of  spore  to  repair  itself  even  if  both  strands  of  the 
DNA  baekbone  are  damaged. 

Sutherland  [97]  proposed  a  stoehastie  model  that  eonsiders  the  probability  that  a  lethal 
damage  ean  be  repaired  to  prediet  a  probability  of  survival  eurve.  Lethal  damage  is  defined 
as  all  lesions  that  ean  eontribute  to  the  inaetivation  of  a  eell  despite  their  ability  to  be 
repaired  or  not.  The  probability  of  a  lethal  damage  oeeurrenee  is  assumed  to  be  given  by 
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the  Poisson  probability  function  and  averaged  repair  probabilities  and  damage  frequencies 
allow  the  calculation  of  an  analytical  solution  in  terms  of  Euler  gamma  functions.  The 
model  also  assumes  the  cell  has  a  limited  repair  capacity.  The  next  section  presents 
the  methodology  utilized  in  this  research  to  develop  a  probability  of  survival  curve  for 
a  population  of  spores. 


5.2  Methodology 

This  research  will  demonstrate  a  stochastic  model  based  on  sampling  from  each  of 
the  spore  population’s  enzyme  ‘fitness’  distributions  in  order  to  determine  the  population’s 
probability  of  survival.  Recall  from  Section  3.8.3  that  each  of  the  six  enzymes  are  assumed 
to  be  represented  by  a  probability  density  function  (PDF)  that  reflects  the  spore  population’s 
initial  enzyme  ‘fitness’  and  that  each  enzyme’s  ‘fitness’  is  mutually  independent.  That  is, 
the  ‘fitness’  associated  with  any  one  enzyme  has  not  effect  or  bearing  on  the  ‘fitness’  of  any 
other  enzyme.  In  addition,  the  initial  ‘fitness’  of  each  enzyme  is  degraded  due  to  oxidative 
damage  caused  by  ROS  within  the  spore  core. 

The  first  step  of  constructing  the  population’s  probability  of  survival  is  sampling  from 
the  cumulative  distribution  function  (CDF)  of  each  damaged  enzyme.  Sampling  from  the 
CDF  allows  the  influence  of  occurrence  of  each  enzyme  to  be  reflected.  The  CDF,  denoted 
by  Ox(.r),  of  the  Uniform  distribution  is  found  by  integrating  Equation  3.34  with  respect  to 
.r.  Ox(.r)  is  given  by 


0,  X  <  a 

fi>x(.^)  =  {  a<x<b 

1,  x>b 

where  Ox(.v)  is  bounded  between  a  and  b.  Figure  5.2  illustrates  the  CDF  of  the  damaged 
UvrB  enzyme  and  demonstrates  values  of  two  samples  from  the  CDF.  It  assumes  a  single 
dose  of  radiation  produced  10  percent  damage  to  the  enzymes.  The  PDF  of  the  UvrB 
enzyme’s  functional  concentration  after  radiation  exposure  has  a  mean  of  336.2  nM  which 
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CDF  of  UvrB  Concentration  (10%  initial  damage)  (nM) 


Figure  5.2:  CDF  and  Sampling  for  UvrB 


is  90  percent  of  the  average  UvrB  concentration  prior  to  irradiation,  a  and  b  are  equal  to 
302.6  nM  and  369.9  nM  respectively.  A  total  of  ten  evenly  spaced  samples  are  taken  from 
the  CDFs  of  each  of  the  six  repair  enzymes.  The  samples  from  the  damaged  enzymes’ 
CDFs  are  used  as  the  initial  enzyme  concentrations  for  the  repair  model  presented  in 
Chapter  4.  Ten  samples  from  each  of  the  six  enzyme  distributions  result  in  a  total  of  one 
million  combinations  of  initial  enzyme  concentrations  that  are  considered  to  determine 
the  concentration  of  viable  DNA,  [Dg]  (7) ,  [Dg]o  (7)),  for  a  given  repair  time,  t.  The 
initial  concentration  of  viable  DNA,  and  the  state  of  the  enzyme  population, 

is  a  function  of  dose,  d.  Dose,  d,  is  equal  to  rC  where  r  is  the  radiation  rate  and 
4  is  the  radiation  exposure  time.  Figure  5.3  portrays  the  estimated  percentage  of  viable 
DNA  (i.e.  the  vertical  axis)  as  a  function  of  time  given  that  ten  percent  of  the  initial  DNA 
concentration  was  damaged.  The  horizontal  axis  represents  time  in  hours.  The  dashed 
red  line  demonstrates  the  DNA  repair  utilizing  the  mean  values  of  the  enzyme  CDFs  as 
the  initial  enzyme  concentrations  and  the  NER  process  with  the  UvrA2Bi  complex.  This 
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output  corresponds  to  the  results  of  the  NER  meehanism  found  in  Seetion  4.6.5  and  Figure 
4.4.  The  blaek  line  represents  repair  over  time  with  the  maximum  eoneentration  values 
sampled  from  eaeh  of  the  enzyme  CDFs  while  the  blue  line  gives  repaired  DNA  over  time 
with  the  minimum  enzyme  eoneentration  values. 


Figure  5.3:  Estimated  Minimum,  Mean,  Maximum  Representations  of  Repaired  DNA 
(10%  initial  damage) 


A  probability  of  kill,  P^,  eurve  for  the  spore’s  population  will  be  eonstrueted  at  three 
different  times  during  the  repair  proeess.  Note  that  probability  of  survival,  Ps  is  equivalent 
to  1  -  Pk-  The  three  separate  times  demonstrate  the  differenee  between  a  population’s 
ehanee  of  survival  when  only  a  short  time  period  (i.e.  50  min)  is  permitted  for  repair  versus 
when  the  repair  proeess  is  allowed  to  proeeed  uninterrupted  for  a  longer  period  of  time  (i.e. 
80  min).  In  addition,  the  probability  of  survival  given  one  hour  of  repair  time  will  also 
be  ealeulated  to  exhibit  the  evolution  of  a  population’s  probability  of  survival  between  the 
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short  and  long  repair  times.  In  order  to  develop  the  kill  curve,  the  concentration  of  viable 
DNA,  is  captured  for  all  of  the  sample  runs  where  tj  represents 

the  different  repair  times  (j  =  1, 2, 3).  The  three  vertical  lines  in  Figure  5.3  depict  the  three 
times  at  which  a  probability  of  kill  curve  will  be  calculated.  Recall  that  at  each  specified 
time,  one  million  data  points  for  [Dg]  (tj,  ^  (7) ,  [T)g]o  (7))  are  captured. 

Histograms  are  constructed  at  each  of  the  three  specified  times  by  dividing  the 
horizontal  axis,  percentage  of  viable  DNA,  into  equal  length  bins  and  drawing  a  rectangle 
with  an  area  proportional  to  the  number  of  observations  in  that  bin.  The  bins  are  labeled 
with  the  minimum  viable  DNA  percentage  from  that  bin.  Observe  from  Figure  5.3  that 
after  50  minutes  of  repair  time,  the  amount  of  viable  DNA  is  roughly  symmetric  about 
the  mean.  This  is  further  demonstrated  in  the  corresponding  histogram  (see  Figure  5.4(a)). 
Histograms  for  the  concentrations  of  viable  DNA  after  one  hour  and  80  minutes  of  repair 
are  also  displayed  in  Figure  5.4.  Note  the  range  of  the  horizontal  axis  of  each  histogram 
is  not  the  same.  Figures  5.4(b)  and  5.4(c)  show  that  with  additional  repair  time,  the 
distribution  of  viable  DNA  becomes  more  skewed  towards  100  percent  repair  of  the  DNA. 

An  estimate  of  the  probabililty  of  kill  is  determined  from  the  distribution  of 
[Dg]  ,  [Dg]o  ^7^^  at  each  tj  by  calculating  its  empirical  distribution  function,  F„(jc), 

defined  by 

Fn{x)  =  -{#Xi<x). 
n 

Here,  n  is  the  number  of  samples  such  that  1  <  z  <  n  and  #A,  is  the  count  of  observations 
of  the  random  variable,  X,  less  than  the  value,  x.  For  this  model,  the  i*  observation  of  the 
random  variable,  A„  is  given  by  [Dg]  (7) ,  [T)g]o  (7)).  and  n  is  one  million.  For  ease  of 
notation,  the  brackets  are  dropped.  The  estimated  cumulative  probability  of  kill,  P^,  for  a 
spore  population  is  approximated  by  its  empirical  CDF  for  n  sufficiently  large. 
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(a)  Histogram  at  50  Min 


(b)  Histogram  at  1  Hr 


Viable  DNA(%) 


(c)  Histogram  at  80  Min 

Figure  5.4:  Histograms  of  Viable  DNA  (%)  at  50  Min,  1  Hr,  and  80  Min  (10%  initial 
damage) 


“While  there  can  be  repair  of  DNA  damage  ...  if  too  much  damage  has  been 
accumulated  during  spore  dormancy,  this  damage  can  overwhelm  the  capacity  of  repair 
systems  and  lead”  to  spore  death  [91].  Let  Dc  denote  the  critical  threshold  such  that  if 
Dg  >  Dc  then  the  spore’s  DNA  has  been  repaired  to  a  viable  level  at  which  it  can  produce  a 
cell  during  germination.  Thus  if  Dg  >  Dc,  the  spore  population  will  survive  and  if  Dg  <  Dc 
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the  population  will  be  die,  i.e. 


K 


Dg  t,  e  ,  <  Z), 


Fn(D,)=- 

n 


#Dg  t,e  -  ,  Dgg  -  <  D, 


Survival  models  found  in  the  literature,  such  as  the  linear-quadratic  model,  relate 
radiation  dose  to  survival.  To  determine  the  explicit  relationship  between  survival  and 
dose,  St(d),  note  that  the  probability  of  survival,  Ps,  is  equivalent  to  1  -  Pk,  therefore 


S,(d)  =  Ps  >  D,|.  (5.1) 

Define  r(J)  to  be  the  support  of  0  {e  the  joint  density  function  given  by  Equation  3.36, 
i.e. 

T(d)  =  {e  :  et  e  [o.9/,(t);U,.^ ,  l.l/.(t);U,,J  ,  /  =  1, . . . ,  b) .  (5.2) 


Let  Q.t(d)  be  the  region  given  by  the  viable  DNA  which  has  reached  the  critical  DNA 
threshold  after  repair  time,  t,  or 

Q,id)  =  :  Dgl^t,e,Dg^  >  D,| .  (5.3) 

Survival  of  the  spore  population,  S  t(d),  is  represented  by  the  intersection  of  r(J)  and  Qf(^): 

St{d)=  J  0(e)Dg|t,e,Dg„|^jjde. 

r(d)nn,(d) 

For  this  research,  the  joint  density  function  is  uniform  for  e  therefore. 


The  results  from  the  probability  of  kill  model  are  portrayed  in  the  next  section. 


5.3  Results 

The  probability  of  kill  of  a  spore  population  is  first  examined  for  the  case  of  ten  percent 
damage  to  its  DNA  and  repair  enzymes  and  utilizes  NER  repair  with  the  UvrA2Bi  complex. 
Assume  the  critical  threshold,  D^,  is  98  percent.  Therefore,  DNA  must  be  repaired  to  98 


90 


percent  of  its  original  capacity  for  the  spore  to  accomplish  necessary  biological  functions. 
This  critical  threshold  is  indicated  by  the  vertical  yellow  lines  in  the  histograms  of  Figure 
5.4.  Notice  that  after  50  minutes  of  repair  time,  the  majority  of  the  population  has  not 
met  this  critical  threshold  (see  Figure  5.4(a)).  The  corresponding  estimated  probability 
of  kill  curve  for  this  spore  population  is  given  by  the  blue  line  in  Figure  5.5.  Again,  the 
vertical  yellow  line  represents  Dc  and  the  horizontal  axis  is  percentage  of  viable  DNA.  The 
probability  of  kill  is  given  by  the  vertical  axis.  Therefore,  given  a  repair  time  of  50  minutes, 
about  89  percent  of  the  spore  population  will  fail  to  meet  the  critical  threshold  of  98  percent 
DNA  viability  i.e. 

50  min,e|-j,90% 

However,  as  shown  by  the  red  line  in  Figure  5.5,  after  one  hour  of  repair,  approximately 
eight  percent  of  the  spore  population  has  failed  to  meet  the  critical  DNA  threshold.  This  is 
also  demonstrated  by  the  histogram  in  Figure  5.4(b).  Also  note  that  80  minutes  of  repair 
time  allows  the  estimated  probability  of  kill  for  a  spore  population  to  be  zero  for  the  critical 
threshold  of  98  percent  viable  DNA.  Thus  the  population  has  an  estimated  100  percent 
probability  of  survival  (demonstrated  by  the  black  kill  curve  in  Figure  5.5). 

Next,  the  effect  of  simultaneous  hydrolysis  damage  on  the  population’s  probability  of 
kill  is  explored.  Figure  5.6  depicts  the  estimated  kill  curves  for  the  same  repair  times  as  the 
previous  graphs  (Figure  5.5)  but  includes  the  possibility  of  continual  hydrolysis  damage  to 
the  DNA  given  the  spores  in  the  population  have  a  5  percent  moisture  content.  Observe 
that  a  50  minute  repair  time  results  in  an  estimated  100  percent  probability  of  kill  for  the 
population  if  the  critical  threshold  of  viable  DNA  is  98  percent.  Even  with  an  hour  of  repair 
time,  about  65  percent  of  the  population  will  fail  to  meet  the  critical  threshold  for  viable 
DNA.  However,  80  minutes  of  repair  time  gives  an  estimated  100  percent  probability  of 
population  survival. 


<  98%  ^  =  89%. 
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Figure  5.5:  Estimated  Probability  of  Kill  Assuming  Repair  with  UvrA2Bi  (10%  initial 
damage) 


Figure  5.6:  Estimated  Probability  of  Kill  Assuming  Repair  with  UvrA2Bi  Including 
Hydrolysis  Damage  with  5%  Moisture  Content  (10%  initial  damage) 
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Figure  5.7  represents  the  probability  of  kill  of  a  spore  population  that  experienees 
eontinual  hydrolysis  damage  with  a  nine  pereent  initial  moisture  eontent.  These  results 
indieate  that  the  estimated  probability  of  kill  of  a  population  that  experieneed  10  pereent 
damage  to  its  DNA  and  repair  enzymes,  is  100  pereent  for  repair  times  less  than  or  equal 
to  an  hour.  Again,  this  is  assuming  a  eritieal  DNA  threshold  of  98  pereent.  Even  with 
80  minutes  of  allowed  repair  time,  there  is  an  estimated  98  pereent  probability  that  the 
population  will  fail  to  meet  the  eritieal  DNA  threshold. 


Figure  5.7:  Estimated  Probability  of  Kill  Assuming  Repair  with  UvrA2Bi  Ineluding 
Hydrolysis  Damage  with  9%  Moisture  Content  (10%  initial  damage) 


In  order  to  examine  the  population’s  survival  dependenee  on  radiation  exposure,  reeall 
from  Seetion  3.8.2,  Equation  3.30,  that  the  eoneentration  of  DNA,  [D](t),  is  a  funetion  of 
the  initial  eoneentration  of  DNA,  [DJr,  and  the  fraetion  of  DNA  remaining  after  some 
radiation  exposure  time,  fnit),  i-e. 

[D](0  =  {DhfD{t).  (5.4) 
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The  dose  dependenee  of  this  model  is  represented  by  exposure  time  whieh  is  related  to 
dose,  d,  by  the  expression  d  =  rte  where  r  is  the  radiation  rate.  Thus  Equation  5.4  beeomes 

[D]|^j  =  [D]r/o|^j.  (5.5) 

As  dose  inereases,  the  fraetion  of  viable  DNA  deereases.  The  enzyme  aetivity  level  as  a 
funetion  of  exposure  time  is  modeled  similarly  to  Equation  5.5  (see  Seetion  3.8.3).  Eigure 
5.8  depiets  the  dose  dependenee  on  repaired  DNA  assuming  no  simultaneous  hydrolysis 
damage  oeeurs.  The  vertieal  axis  represents  the  pereentage  of  viable  DNA  while  the 
horizontal  axis  gives  time  in  hours.  If  both  enzymes  and  DNA  reeeive  a  radiation  dose 
that  damages  their  initial  eoneentrations  by  20  pereent,  the  neeessary  DNA  repair  time  will 
be  inereased  as  shown  the  by  the  red  line  in  Eigure  5.8.  Utilizing  the  mean  value  from  eaeh 
of  the  six  enzyme  ‘fitness’  CDEs  that  were  damaged  by  20  pereent  results  in  an  estimated 
repair  time  of  approximately  3  hours.  This  is  more  than  twiee  as  long  as  the  repair  of  DNA 
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which  only  experienced  10  pereent  initial  damage  to  its  DNA  and  the  state  of  its  enzymes 
(shown  by  the  blue  line).  The  blaek  line  in  Figure  5.8  represents  repair  of  DNA  given  that 
30  pereent  of  its  initial  DNA  and  repair  enzymes  were  damaged.  After  four  hours  of  repair, 
the  DNA  still  has  not  been  restored  to  95  pereent  of  its  original  state. 

Next,  the  probability  of  kill  for  a  population  that  experieneed  a  radiation  dose  that 
damaged  its  DNA  and  the  state  of  its  enzymes  by  20  pereent  will  be  determined.  Figure  5.9 
portrays  the  estimated  pereentage  of  repaired  DNA  (i.e.  the  vertieal  axis)  as  a  funetion  of 
time  given  that  20  pereent  of  the  initial  DNA  eoneentration  was  damaged  and  there  exists  no 
simultaneous  hydrolysis  damage.  The  horizontal  axis  represents  time  in  hours.  The  dashed 
red  line  demonstrates  the  DNA  repair  utilizing  the  mean  values  of  the  enzyme  CDFs  as 
the  initial  enzyme  eoneentrations  and  the  NER  proeess  with  the  UvrA2Bi  eomplex.  The 
blaek  line  represents  repair  as  a  funetion  of  time  using  the  maximum  eoneentration  values 


Figure  5.9:  Estimated  Minimum,  Mean,  Maximum  Representations  of  Repaired  DNA 
(20%  initial  damage) 
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sampled  from  each  of  the  enzyme  CDFs  while  the  blue  line  gives  repaired  DNA  over  time 
with  the  minimum  enzyme  concentration  values.  For  all  three  cases,  the  DNA  repair  time 
required  is  more  than  twice  as  long  as  the  DNA  repair  time  when  only  10  percent  of  the 
DNA  and  enzymes  were  damaged. 

Figure  5. 10  illustrates  the  effect  of  an  increased  dose  on  the  population’s  probability  of 
kill.  These  kill  curves  are  calculated  for  a  population  that  experienced  20  percent  damage 
to  its  DNA  and  repair  enzymes  but  did  not  encounter  any  simultaneous  hydrolysis  damage. 
Note  that  these  kill  curves  represent  different  repair  times  than  the  previous  probability  of 
kill  curves  (Figure  5.5).  With  the  increased  damage  to  DNA  and  enzymes,  an  estimated 
100  percent  of  the  population  will  fail  to  meet  the  critical  threshold  of  98  percent  viable 
DNA  after  80  minutes  of  repair  time  (refer  to  the  blue  line  in  Figure  5.10).  A  population’s 
estimated  probability  of  kill  is  still  99  percent  after  two  hours  of  repair  as  shown  by  the  red 
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Viable  DNA(%) 

Figure  5.10:  Repair  with  UvrA2Bi  (20%  initial  damage) 
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line  in  Figure  5.10.  The  black  line  indicates  that  approximately  1 8  percent  of  the  population 
will  fail  to  meet  the  critical  threshold  of  98  percent  viable  DNA  after  two  and  a  half  hours 
of  repair. 

5.4  Summary 

This  chapter  combined  damage  to  a  spore’s  DNA  and  enzymes  with  DNA  repair 
to  develop  and  explore  a  probability  of  kill  model  for  a  population  of  Bacillus  spores. 
A  critical  threshold  of  DNA  viability  was  assumed  and  the  effect  of  repair  time  on  a 
population’s  probability  of  survival  was  examined.  An  increase  in  allowable  repair  time 
decreases  the  population’s  kill  probability.  However,  simultaneous  hydrolysis  damage 
increases  the  probability  of  kill  and  a  spore  population’s  chance  of  surviving  decreases  as 
the  moisture  content,  and  thus  the  hydrolysis  damage,  increases.  In  addition,  a  larger  dose 
of  radiation  increases  the  damage  to  both  DNA  and  repair  enzymes  which  significantly 
decreases  a  spore  population’s  survivability  even  with  increased  repair  time.  Finally,  an 
expression  for  population  survival  as  a  function  of  radiation  dose  was  found. 
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VI.  Conclusion 


This  research  models  the  inactivation  of  a  population  of  Bacillus  spores  after 
irradiation  by  coupling  radiation  dose,  damage,  and  repair  in  order  to  characterize  the 
population’s  probability  of  survival.  Indirect  damage  to  a  Bacillus  spore’s  DNA  and 
essential  enzymes  was  modeled  via  an  ionizing  radiation  exposure.  The  irradiation  of  water 
produces  radicals  which  react  with  biological  material  as  they  diffused  into  the  spore.  As  a 
result,  these  reactive  oxygen  species  (ROS)  damaged  the  spore’s  DNA  and  repair  enzymes 
within  the  core.  The  extent  of  this  damage  hampered  the  ability  of  the  spore’s  DNA  to  be 
repaired  and  thus  the  spore’s  germination  and  outgrowth  capability.  All  of  these  factors 
effect  a  population’s  probability  of  survival. 

The  production  and  diffusion  of  ROS  and  their  subsequent  reactions  within  the  spore’s 
core  were  shown  to  be  dependent  on  the  radiation  exposure.  This  research  considered 
a  20  Mev,  1  ns  pulse  of  radiolysis  impinging  on  an  inert  solid  plate  having  a  thin  layer 
of  adsorbed  water  and  an  adsorbed  spore;  simulating  inactivation  of  spores  adsorbed  on 
surfaces  in  a  humid  environment.  This  produced  a  radical  yield  on  the  outer  boundary  of 
the  spore  which  was  in  an  ionizing  particle’s  spur.  The  reaction-diffusion  of  these  ROS  into 
the  spore  core  resulted  in  damage  to  the  spore’s  DNA  such  as  DNA  base  modifications  and 
abstraction  of  molecules  from  the  DNA  strand.  The  average  concentration  of  reactants  in 
the  core  led  to  approximately  10  percent  of  the  spore’s  DNA  experiencing  damage  due  to 
the  irradiation  considered  in  this  research.  The  spore’s  repair  enzymes  were  also  degraded 
due  to  ROS  reactions. 

The  repair  mechanisms  of  Bacillus  spores  must  be  considered  prior  to  examining  a 
spore  population’s  ability  to  survive.  This  research  developed  a  model  for  nucleotide 
excision  repair  (NER)  which  is  the  primary  mechanism  by  which  indirect  damage  via 
ionizing  radiation  is  repaired.  The  NER  process  repairs  a  section  of  damaged  DNA  and 
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it  is  responsible  for  the  removal  of  a  broad  group  of  strueturally  unrelated  DNA  lesions  to 
inelude  SSBs  and  the  interaetions  of  damaged  bases.  Variations  of  the  NER  pathway  were 
eonsidered  as  well  as  an  extension  that  ineluded  the  possibility  of  simultaneous  hydrolysis 
damage  during  the  repair  proeess.  Beeause  hydrolysis  damage  oeeurs  simultaneously  with 
repair,  the  level  of  viable  DNA  may  never  reaeh  the  neeessary  threshold  for  outgrowth.  In 
addition,  the  length  of  required  repair  time  and  the  rate  of  repair  is  dependent  upon  the 
moisture  eontent  of  the  spore. 

The  probability  of  kill  for  a  spore  population  was  examined  under  the  assumption 
there  exists  some  eritieal  threshold  of  viability  that  the  DNA  must  be  repaired  to  in  order 
for  the  population  to  eonduet  neeessary  biologieal  proeesses  for  survival.  In  addition,  a 
eertain  initial  enzyme  ‘fitness’  of  the  population  was  taken  into  aeeount.  As  expeeted, 
the  population’s  probability  of  kill  deereases  as  allowed  repair  time  inereases.  However, 
simultaneous  hydrolysis  damage  has  a  signifieant  effeet  on  a  population’s  probability  of 
kill  and  inereasing  moisture  eontent  of  the  population  further  inereases  its  kill  probability. 
Also,  inereasing  the  radiation  dose  inereases  the  damage  via  ROS  to  both  the  spore’s  DNA 
and  enzymes.  Subsequently,  the  population’s  probability  of  kill  is  inereased  when  radiation 
dose  inereases.  This  researeh  provides  a  eharaeterization  of  a  spore’s  probability  of  kill 
depending  on  radiation  exposure,  ehemieal  reaetions  within  the  spore,  and  repair  time. 

6.1  Contributions 

Contributions  to  the  field  of  work  provided  by  this  researeh  inelude  a  prokarytie  NER 
model.  To  date,  only  the  prokaryotie  base  exeision  repair  (BER)  and  non-homologous 
end-joining  proeesses  have  been  mathematieally  modeled  and  even  these  models  utilized 
eukaryotie  enzyme  rate  kineties.  In  addition,  this  researeh  provided  the  first  deseription 
of  radiation  damage  in  a  humid  environment  in  whieh  hydrolysis  eontributes  to  the 
inaetivation  of  the  spore.  Einally,  a  proeess  by  whieh  dose  is  direetly  related  to  probability 
of  survival  was  presented.  This  was  done  by  eoupling  radiation  dose,  damage  to  DNA 


99 


and  enzymes,  and  DNA  repair  while  assuming  a  certain  population  initial  enzyme  ‘fitness’ 
probability  density  function. 

6.2  Future  Work 

During  development  of  the  damage,  repair,  and  probability  of  kill  models  contained 
in  this  research,  the  following  areas  were  identified  as  requiring  further  analysis.  The 
damage  model  assumed  that  there  exist  no  significant  permeability  barriers  within  the 
spore  that  would  prevent  the  movement  of  reactants  from  the  outer  boundary  to  the  core. 
As  discussed  in  Section  2.1.1  the  exosporium,  the  spore  coat,  the  outer  membrane,  the 
cortex,  and  the  germ  cell  wall  are  not  considered  to  be  permeability  barriers  against  small 
molecules.  However,  the  inner  membrane  of  the  dormant  spore  might  prevent  a  portion  of 
the  reactants  from  reaching  the  core  and  thus  the  spore’s  DNA.  Therefore  the  diffusion  of 
reactants  through  this  membrane  could  be  incorporated  into  the  damage  model. 

The  radiation  damage  model  considered  two  classes  of  reactants.  Ri  represents 
the  class  that  reacts  quickly  with  the  O2  present  in  the  adsorbed  water  and  therefore 
instantaneously  transforms  into  other  species.  e~  is  an  example  of  this  class  of  reactants. 
R2  represents  the  class  that  does  not  react  as  fast  as  the  Ri  class  and  therefore  does  not 
get  transformed  as  quickly.  This  allows  the  class  of  R2  reactants  to  diffuse  to  the  outer 
boundary  of  the  spore.  OH  belongs  to  this  class  of  reactants.  However,  this  model  may  be 
oversimplified.  After  examining  the  rate  coefficients  in  Table  3.2  of  OH  with  other  species, 
a  third  class  of  reactants  should  be  considered.  Specifically,  the  reaction  of  OH  with  O^ 
has  a  reaction  rate  which  is  almost  equivalent  to  the  diffusion  controlled  rate.  However,  the 
reaction  of  OH  with  H2O2  has  a  reaction  rate  which  is  much  slower  (3.8  x  10^ 
Therefore,  a  damage  model  that  breaks  up  the  current  definition  of  class  2  reactants  into  two 
separate  categories  thereby  defining  three  different  classes  of  reactants  could  be  considered. 
The  current  R2  class  would  be  divided  into  a  category  of  reactants  that  will  reach  the  outer 
boundary  of  the  spore  and  a  category  of  reactants  that  would  react  with  other  species  as 
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it  diffuses  through  the  water  layer  between  the  souree  and  the  spore.  Furthermore,  there 
exists  a  possible  need  for  four  reaetant  elasses  in  order  to  model  the  reaetion-diffusion  of 
the  OH,  H,  e~,  and  H2O2  speeies. 

In  addition,  the  eoneentration  of  reaetants  was  assumed  to  be  uniformly  impressed 
on  the  outer  boundary  of  the  spore.  However,  if  the  substrate  is  irradiated  at  only  one 
loeation  as  in  Figure  3.4  and  the  spore  is  loeated  relatively  close  to  the  radiation  source,  the 
reactants  would  not  be  uniformly  distributed  within  the  substrate.  A  greater  concentration 
of  reactants  would  be  expected  on  one  side  of  the  spore  compared  to  the  other  side  of  the 
spore.  For  the  case  considered  in  this  research,  the  radiation  events  would  be  randomly 
distributed  around  the  spore.  Integrating  over  this  distribution  of  event  locations  could  be 
an  approach  to  determine  the  distribution  of  reactants  on  the  outer  boundary  of  the  spore. 

Certain  values  given  in  Section  4.6.5  for  reaction  rate  parameters  of  the  NER  are  not 
well  known.  For  example,  the  rate  at  which  polymerase  I  lays  down  corrected  bases  found 
in  the  literature  applies  to  BER.  During  BER,  only  one  base  is  replaced  while  NER  requires 
multiple  bases  to  be  replaced.  Sensitivity  analysis  on  all  of  the  repair  rate  parameters  should 
be  completed.  Also,  the  NER  model  and  thus  the  probability  of  survival  model  presented 
in  this  research  assumed  that  no  repair  molecules  were  produced  after  the  initiation  of 
germination.  Euture  work  should  examine  the  possibility  of  new  repair  enzyme  creation 
from  undamaged  or  repaired  DNA. 

The  initial  enzyme  ‘fitness’  presented  in  Section  3.8.3  was  modeled  via  a  Uniform 
PDE.  This  representation  is  conservative  and  the  effects  of  various  density  functions  should 
be  examined.  In  addition,  the  utilization  of  a  continuous  PDE  was  necessary  due  to  the 
modeling  of  the  repair  process  via  differential  equations.  However,  the  low  concentrations 
of  the  UvrA  and  UvrC  enzymes  might  be  more  accurately  modeled  with  a  discrete 
distribution  or  a  mixture  distribution.  Such  an  approach  would  require  the  NER  process 
to  be  modeled  as  a  discrete  state  difference  equation. 
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Further,  the  histograms  given  in  Figure  5.4  whieh  eharaeterize  the  distribution  of  the 
repaired  DNA  over  time  indieate  that  it  may  be  possible  to  model  the  repaired  DNA  with 
a  PDF  as  a  funetion  of  time.  Reeall  that  the  distribution  of  repaired  DNA  is  related  to 
the  initial  enzyme  ‘fitness’  PDFs.  Modeling  the  repaired  DNA  as  a  PDF  over  time  would 
eliminate  the  need  to  eonduet  million-run  simulations  whieh  require  extensive  eomputing 
power  and  time.  However,  in  order  to  determine  the  density  funetion  of  the  repaired  DNA, 
additional  samples  from  eaeh  of  the  CDFs  of  the  six  repair  enzymes  should  be  eonsidered. 
This  researeh  took  ten  evenly  spaced  samples  from  each  enzyme’s  CDF.  In  order  to  gain 
the  fidelity  needed  to  calculate  the  repaired  DNA  distribution,  one  hundred  samples  from 
each  enzyme  CDF  should  be  taken  into  account. 

Finally,  further  analysis  is  necessary  to  compare  this  research  to  real  data.  Experiments 
conducted  with  live  spores  would  allow  comparison  with  data  used  in  the  mathematical 
models  presented  in  this  research.  This  research  demonstrated  a  method  to  compute 
probability  of  kill  curves  but  many  assumptions  were  made  during  their  development  such 
as  the  critical  DNA  threshold  value.  These  assumptions  should  be  examined  further  in 
order  to  incorporate  uncertainty  and  replicate  with  experimental  results.  Initial  results 
of  this  cursory  look  at  the  survival  curves  is  promising.  Functional  data  analysis  to 
refine  the  enzyme  distributions  should  be  completed  In  addition,  replication  of  the  method 
with  experimental  results  and  comparison  with  other  methods  of  probability  of  survival 
estimation  could  be  examined. 
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Appendix:  Steady  State  Solution  to  the  Mathematical  Repair  Model 


This  appendix  presents  the  steady  state  solution  to  the  NER  mathematical  model 
presented  in  Sections  4.6.1  and  4.6.4.  This  steady  state  solution  is  obtained  when  the 
time  evolution  of  the  system  vanishes.  Therefore,  the  differential  equations  are  set  to  zero 
producing  the  algebraic  system: 

0  =h  ([A2]o  -  [A2]  -  [A2B])  +  k-AAiB]  -  h[A2\m  (A.l) 

0  ([^]o  -  [//])  [P]  +  k-dAiB]  -  kdA2][B]  (A.l) 

0=k5([C]o-[C])[H] 

-  k,  {[fi]o  +  [A2]  -  [A2]o  -  [B]  -  ([C]o  -  [C])  -  ([H]o  -  [//])}  [C]  (A.3) 

0  =ke  ([H]o  -  [//])  [P]  -  h  ([C]o  -  [C])  [H]  (A.4) 

0  =^7  ([P]o  -  m)  [L]  -  ke  ([//]o  -  [H])  [E]  (A.5) 

0  =k,  ([L]o  -  [L])  -  kj  ([E]o  -  [P])  [L]  (A.6) 

0  =  -  k2[A2B][D]  +  k*[Dg]  (A.7) 

0  =ki[A2][B]  -  k_i[A2B]  -  k2[A2B][D]  (A.8) 

0=fc8([^]o-[^])-fc*[^g].  (A.9) 

For  convenience  we  set 

x  =  k*Dg,  y  =  [A2B],  and  z  =  [A2],  (A. 10) 

adopt  the  notation  zo  =  [A2]o,  and  drop  the  brackets.  Upon  substitution  and  some  algebraic 
reduction,  we  obtain  a  new  set  of  equations: 

ks(Lo-L)-x  =  0  (A.ll) 

kj(Po  -  P)L  -  X  =  0  (A.  12) 

h(HQ-H)P-x  =  0  (A.  13) 
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ksiCo  -  C)H  -  X  =  0  (A.  14) 

hizo  -z)-hy-  x  =  0  (A.  15) 

k\zB  -  k-\y  -  X  =  0  (A. 16) 

k2yD-x  =  0  (A.  17) 

k4C[(Bo  -B)-  izo  -z)-  (Co  -C)-{Ho-H)]-x  =  0  (A.18) 


Equation  A.  11  is  from  A.9  and  Equations  A.  12  through  A.  14  are  obtained  by  sequential 
substitution  of  this  result.  Equations  A.l,  A.2,  and  A. 7  are  Equations  A. 15,  A. 16,  and 
A.  17  respectively.  Einally,  Equation  A.18  is  Equation  A.3  with  the  replacement  given  from 
Equation  A.  14.  The  results  from  Equation  A. 8  are  not  included  as  they  only  reproduce 
Equation  A.  16.  An  additional  equation  is  necessary  to  complete  the  algebraic  system.  It  is 
attained  from  the  DNA  conservation  requirement  that 

Dj  =  Dg  +  Z)  +  D(- 

as  described  in  Section  4.6.4  with 

Dc  =  A2BD  +  HBD  +  CBD  +  BD  +  PD  +  LD. 

Using  the  enzyme  conservation  equations  of  Section  4.6.1,  we  obtain 

Dr  =  Z),  +  D  +  (5o  -  5)  +  (Eq  -  P)  +  iU  -  L)  -  {A2B]. 

Substituting  the  definitions  of  Equation  A.  10  results  in 

^  +  Z)  +  (5o  -  5)  +  (Do  -  E)  +  (Lo  -  L)  -  Dr  -y  =  0  (A.19) 

Now  Equations  A.  1 1  through  A.19  form  a  complete  set  of  algebraic  equations  to  determine 
the  steady  state  behavior  of  the  NER  process. 

Begin  by  examining  the  case  where  no  hydrolysis  damage  is  occurring  (i.e.  k*  =  0) 
which  implies  x  =  0.  Then  from  Equations  A.l  1  through  A.  14  it  is  clear  that 

L  =  Lo,  P  =  Po,  H  =  Ho,  and  C  =  Cq.  (A.20) 
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Consequently,  over  time,  the  concentrations  of  UvrC,  helicase  II,  polymerase  I,  and  ligase 
will  approach  their  initial  concentrations.  From  Equation  A.  17,  we  find  either  y  =  0  or 
D  =  0.  If  y  =  0  then  Equation  A.  15  implies  z  =  zo-  Substituting  this  along  with  Equation 
A.20  into  Equation  A.  18  leads  to  B  =  Bq.  Now  Equation  A.  16  becomes  zB  =  zqBq  4^  0 
which  is  a  contradiction.  Therefore  y  4^  and 


D  =  0  (A.21) 

which  implies  complete  DNA  repair  has  taken  place.  Now  Equation  A.  17  is  satisfied.  Erom 
Equation  A.  15  and  subsequently  from  Equation  A.  18  we  have 

z  =  zo  -  y  and  B  =  B^-y.  (A. 22) 

Note  that  Equation  A.  16  is  the  equilibrium  equation  implying  that  the  concentrations  of 
UvrA,  UvrB,  and  the  UvrA2Bi  complex  will  all  exist  in  their  equilibrium  state  given  by 
the  equation  k-i\A2B\  =  ki[A2\\_B\  where  ^  is  the  equilibrium  constant  [85].  Substituting 
Equation  A. 22  into  the  equilibrium  equation  produces 


^i(zo  -  y){BQ  -  y)  -  k-iy  =  0. 

Solving  this  quadratic  equation  for  y  with  d  =  ^  admits  two  solutions, 
y  =  2  +  zo  +  Bq)  +  V(d  +  Zo  +  -  4zoBoj . 

Because  as  Bq  ^  0  requires  [A2B]  =  y  ^  0,  it  follows  the  minus  sign  must  be  chosen. 
Therefore 

[A2B]  =  y  =  1  [(d  +  zo  +  Bo)  -  V(d  +  zo  +  Bo)2-4zoBo]  (A.23) 

Thus  the  steady  state  solution  with  no  hydrolysis  damage  included  is  given  by  Equations 
A.20  through  A.23. 
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Next  explore  the  case  where  k*  0.  Using  Equations  A.  11  through  A.  19  and  solving 


for  L,  P,  B,  H,  C,  y,  and  z  in  terms  of  v:  and  D  gives 
xk^ik-i  +  Dk2) 


B  =  - 


C  =  Co- 
H  =  Ho- 


k\  \_x{k2,  +  Dk2)  -  Dk2k2Zo] 

xke  x{h  +  Poh)  -  LoPokjks 


ks  -x^k-j  +  x(HQk(,k8  +  L^k-ik^  +  HqP okf^k-i)  —  HqLqP ()k(,kjk^ 

xkj  X  —  Lok^ 


ke  x(kg  +  P okj)  -  LqP okjkg 


L  =  Lo  -  - 

P  =  Po  + 


X 

h 

X 


kj  X  —  L{)k^ 


y  = 


Dk2 

Z  =  Zo-  X 


1  1 

Dk2  kj, 


(A.24) 

(A.25) 

(A.26) 

(A.27) 

(A.28) 

(A.29) 

(A.30) 


Notice  that  C,  H,  L  and  P  depend  only  on  x  while  B,  y,  and  z  depend  on  both  x  and  D. 
These  results  are  substituted  into  Equations  A.  18  and  A.  19  to  obtain  a  system  of  equations 
for  v:  and  D. 


^4  C 


xk^ 


x(ks  +  P  okj)  —  LqP  okjks, 


ks  —x^kj  +  x(Hok(,ks.  +  Lokjk^  +  HqP  ok^kj)  —  HqLqP  ok^kjks, 


Bn  + 


xksik-i  +  Dk2) 


-  X 


1 


^  1  \  / xkj 


X  -  Lnh 


o^-s 


ki  {x{ks  +  Dk2)  -  Dk2hzo\ )  \Dk2  ksj  \  ke  x{h  +  Poh)  -  LnPnkih 
xke  x(ks  +  Pokj)  -  LoPokjk^ 


-x^kj  +  x(Hnk(,k^  +  Lnkjk^  +  HqP ok^kj)  -  HqLqP ok(,kjks 


X  =  0 


(A.31) 


jc  /d  ,  xk2{k-i+Dk2) 

-  +  Zy  +  I  1JC\  "T  - 

k*  \  k\  {x{k2  +  Dk2)  -  Dk2k2,ZQ\ 


kj  X  —  Lnk^ 


+  I  —  1  —  Dt  — 


Dko 


(A.32) 


Eurther  progress  with  these  equations  must  be  made  numerically. 
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